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ABSTRACT 


The focus of present research work is on the modeling of the event-based rainfall- 
runoff (R-R) process using a variety of conventional and new techniques. Conceptual 
rainfall-runoff (CRR) models were used as conventional techniques and Artificial 
Neural Network (ANN) was used as a new techniques. Among the CRR model four 
widespread conceptual models, the Unit hydrograph model, Nash conceptual model, 
Clark conceptual model and simple non-linear Tank model, were developed. Among 
the ANNs, twenty different architectures having single hidden layer were explored. 
An attempt has also been made to develop new guidelines and methodologies for the 
development of the ANN models. The back propagation algorithm of MATLAB 
toolbox was used to develop ANN models in this study. The stream flow and rainfall 
data from Kentucky River Basin (KRB), USA were employed to develop all these 
models. The worth of multi-criteria validation for evaluating model consistency was 
emphasized. All models were capable to simulate runoff adequately after calibration. 

Among the CRR models, UH model was found to perform best during both 
calibration and validation. ANN 3-10-1 was found best the ANN model when the data 
from many storms are combined into one training set. The evaluation results in terms 
of standard statistical parameters between CRR model simulation result and ANN 
models were found to be comparable. It was observed that rather than relying on a 
few performance statistics such as correlation coefficient, efficiency and root mean 
square error, a wide variety of error statistics including average absolute relative error, 
normalized mean bias error, and threshold statistics should be used to evaluate the 
predictive capability of the developed models. The use of such methods enables 
assessment of the reliability of model predictions. It also supports the further 
development of models by identification of weak parts and evaluation of 
improvements. 


Ill 



ACKNOWLEDGEMENT 


I take this opportunity to express my deep sense of gratitude to Dr. Ashu Jain 
for his constant inspiration, support and encouragement given to me from time to 
time. His guidance, cooperation and meticulous scientific attitude inspired and 
immensely benefited me for the completion of this work. It was indeed a great 
privilege for me to have worked under his supervision. I also pay my sincere gratitude 
to Dr. T. Gangadharaiya, Dr. Rajesh Srivastva, Dr. Bithin Datta and Dr. Pranab 
Mahapatra for their encouragement and affection throughout my stay in IIT Kanpur. 

I am fortunate to have nice and so helpful friends at IITK. In fact it was their 
support, inspiration and love that made my stay memorable and delightful here. Their 
concern, love and support during crisis ai'e some of the things that I will never forget 
in my life. To list them is a Herculean task for me. Simply words are not sufficient to 
acknowledge their kind help. I wish to extend my thanks to all the friends in 
Hydraulics and water resources engineering group in IIT-Kanpur. I acknowledge the 
help received from my seniors, juniors and my batch mates. I enjoyed partying, 
joking, discussing problem with them. In their company, I always relaxed from the 
tension of work and enjoyed most. I would like to thank all of the HWRE lab 
members. 

I would like to express my deepest gratitude to my parents for their love, care 
blessings, throughout and unconditional in my whole academic career. Thanks to all 
people who have helped me directly or indirectly throughout the thesis work and have 
helped me in completing the job successfully. 

Last but not least, I would like to express my thanks and gratitude to the 
almighty, the most beneficent and merciful, for granting me an opportunity to be here 
in world’s one of the best educational institute. I will cherish its memories lifelong. 

July 15, 2004 Subhash Chandra 

Indian Institute of Technology 
Kanpur 


IV 



CONTENTS 


ABSTRACT (iii) 

LIST OF FIGURES (vii) 

LIST OF TABLES (X) 

1. Introduction 

1.1 General 1 

1 .2 Statement of Problem 1 

1.3 Research Objectives 4 

i .4 Organization of Thesis 4 

2. Literature Review 

2.1 General 6 

2.2 CRR Models 6 

2.3 ANN Models 9 

3. Artificial Neural Network 

3.1 General 12 

3.2 The Artificial Neuron 12 

3.3 Architecture of ANNs 13 

3.4 Learning 14 

3.4. 1 The Back-propagation •. 14 

3.5 The ANN Tool Box 16 

4. Mlodel Development 

4.1 General 21 

4.2 Study Area And Data 21 

4.3 Model Development 26 

4.3.1 Unit Hydrograph Model 28 

4. 3. 1.1 Calculation of (!> - Index 30 

4.3.2 Nash Conceptual Model 35 

4.3.3 Clark Model 41 

4.3.4 Nonlinear Tank Model 45 

V 



4.3.4. 1 Determination of Parameters 46 

4. 3.4. 2 Calculation ofExcess Rainfall (R) 47 

4.3.5 ANN Model 50 

4.3.5. 1 Selection of Input and Outputs 50 

4. 3. 5. 2 Training of ANN Models 51 

5. Results and Discussions 

5.1 General 62 

5.2 Statistical Parameters 62 

5.3 Discussion of Results 66 

5.3.1 Results for Combined Input 66 

5.3.2 Results for Category-wise Input 67 

5.3.3 Analysis of Various Performance Statistics 70 

6. Summary, Conclusions and Scope of Future Work 92 

APPENDIX A 

Performance of Various Models in Terms of Standard Statistical Parameters 

APPENDIX B 

Performance of Various Models- Scatter Plots 
REFERENCES 


VI 



LIST OF FIGURES 


Figure 3.1 Model of an Artificial Neuron 13 

Figure 3.2a MATLAB window 16 

Figure 3.2b Neural Network Tool 17 

Figure 3.2c Network type and different function 18 

Figure 3.3a Log sigmoid Transfer function 18 

Figure 3.3b Tan sigmoid Transfer function 19 

Figure 3.3c Linear Transfer function 19 

Figure 3.4a 2-5-1 Architecture 19 

Figure 3.4b Training Window 20 

Figure 3.4c Plot- Error v/s Training cycle 20 

Figure 4. 1 Kentucky River Basin 22 

Figure 4.2 Rainfall and streamflow for the storm of SI (May 1960) 23 

Figure 4.3 Rainfall and streamflow for the storm of S2 (May 1961) 23 

Figure 4.4 Rainfall and streamflow for the storm of SIO (April 1970) 24 

Figure 4.5 Rainfall and streamflow for the storm of S13 (May 1972) 24 

Figure 4.6 Relationship Between Effective Rainfall and Peak Flow from 

Calibration Storms 25 

Figure 4.7 Relationship Between Effective Rainfall and Peak Flow from 

All 13 Storms 25 

Figure 4. 8 Model development procedures 27 

Figure 4.9 Rainfall-Runoff Processes UH Conceptual Model 29 

Figure 4.10 Behaviour of (f) -Index With Ground Water Flow Magnitude 31 

Figure 4. 1 1 Behaviour of -Index With Total Rainfall (2-days previous) 32 

Figure 4. 12 Behaviour of ^ -Index With Total Rainfall (3-days previous) 32 

Figure 4.13 Behaviour of ((> -Index With Total Rainfall (4-days previous) 32 

Figure 4. 14 Behaviour of ^ -Index With Total Rainfall (5-days previous) 32 

Figure 4. 15 Calibrated UHs 34 

Figure 4. 16 NASH Model-Cascade of n Linear Reservoirs 37 

Figure 4.17 UHs for Nash model 39 

i 

Figure 4. 18 Isochrones and TA histogram for a watershed 41 


VII 



Figure 4.19 UHs for Clark model 45 

Figure 4.20 Determination of Parameters of Tank Model: 

kg, ks, ko, Qg and Qs 47 

Figure 4.21 Schematic Plot of Recession Curves on Semi logarithmic Graph 48 

Figure 4.22 Auto-correlation function v/s Lag 51 

Figure 4.21 Plot of Acceptable level of global error and Difference Between 
Training and Testing AARE (training S6) from ANN 2-5-1 Model. ..55 
Figure 4.15 Plot of Acceptable level of global error and Difference Between 

Training and Testing E (training S6) from ANN 2-5-1 Model 56 

Figure 4.16 Difference In Performance During Training And Testing For 

ANN 2-5-1 Model For Different Initial Set Of Weights 56 

Figure 4.17 Difference In Performance During Training And Testing For 

ANN 2-5-1 Model For Different Initial Set Of Weights 57 

Figure 4.27 Scatter Plot of Observed and Difference Between Observed 

and Predicted Streamflow from UH Model During Validation 57 

Figure 4.28 Scatter Plot of Observed and Difference Between Observed 

and Predicted Streamflow from NASH Model During Calibration.... 58 
Figure 4.29 Scatter Plot of Observed and Difference Between Observed 

and Predicted Streamflow from NASH Model During Validation 58 

Figure 4.30 Scatter Plot of Observed and Difference Between Observed and 

Predicted Streamflow from ANN 3-10-1 Model During Training 59 

Figure 4.31 Scatter Plot of Observed and Difference Between Observed and 

Predicted Streamflow from ANN 3-10-1 Model During Testing 59 

Figure 4.32 Behaviour of el With Flow Magnitude from UH Model 60 

Figure 4.33 Behaviour of el With Flow Magnitude from ANN 2-5-1 Model 60 

Figure 4.34 Behaviour of e2 With Flow Magnitude from UH Model 61 

Figure 4.35 Behaviour of e2 With Flow Magnitude from ANN 2-5-1 Model 61 

Figure 5. 1 Scatter Pot for UH-c Model During Calibration 78 

Figure 5.2 Scatter Pot for UH-c Model During Validation 78 

Figure 5.3 Scatter Plot for ANN 2-8-lc During Training 79 

Figure 5.4 Scatter Plot for ANN 2-8-lc During Testing 79 

Figure 5.5 Scatter Plot for ANN 3-10-lc During Training 80 

Figure 5.6 Scatter Plot for ANN 3-10-lc During Testing 80 

Figure 5.7 UH Model Validation from Category-I Storms 81 


VIII 



Figure 5.8 UH Model Validation from Category-II Storms 81 

Figure 5.9 UH Model Validation from Category-Ill Storms 82 

Figure 5.10 UH Model Validation from Category-IV Storms 82 

Figure 5.1 1 Scatter Plot for Category-I storms from ANN 2-5-1 Model 

During Training 83 

Figure 5.12 Scatter Plot for Category-I storms from ANN 2-5-1 Model 

During Testing 83 

Figure 5.13 Scatter Plot for Category-II storms from ANN 2-6-1 Model 

During Training 84 

Figure 5.14 Scatter Plot for Category-II storms from ANN 2-6-1 Model 

During Testing 84 

Figure 5.15 Scatter Plot for Category-III storms from ANN 2-10-1 Model 

During Training 85 

Figure 5.16 Scatter Plot for Category-III storms from ANN 2-10-1 Model 

During Testing 85 

Figure 5.17 Scatter Plot for Category-IV storms from ANN 3-6-1 Model 

During Training 86 

Figure 5.18 Scatter Plot for Category-IV storms from ANN 3-6-1 Model 

During Testing 86 

Figure 5.19 During Calibration/Training Combined Input 87 

Figure 5.20 During V alidation/Testing Combined Input 87 

Figure 5.21 Category-I Storm During Calibration/Training 88 

Figure 5.22 Category-I Storm During Validation/Testing 88 

Figure 5.23 Category-II Storm During Calibration/Training 89 

Figure 5.24 Category-II Storm During Validation/Testing 89 

Figure 5.25 Category-III Storm During Calibration/Training 90 

Figure 5.26 Category-III Storm During Validation/Testing 90 

Figure 5.27 Category-IV Storm During Calibration/Training 91 

Figure 5.28 Category-IV Storm During Validation/Testing 91 


IX 



LIST OF TABLES 


Table 4.1 Runoff depth and (p index for the calibration storms 30 

Table 4.2 Computation of 24-hour UH by UH theory 33 

Table 4.3 Calibrated UHs from UH model 33 

Table 4.4 Runoff Depth And (p Index For Validation Storms 34 

Table 4.5 Calculation of M/, and M,, of Nash’s Model 38 

Table 4.6 Calculation of and M^^of Nash’s Model 38 

Table 4.7 Nash’s Model Parameters- n and K 39 

Table 4.8 Calibrated UHs from NASH Model ...40 

Table 4.9 Calculation 24-hour UH by NASH model 40 

Table 4.10 CLARK Model Parameters- Ci, C 2 and K 43 

Table 4 . 1 1 Calculation of UH Clark’s method 44 

Table 4.12 Calibrated UHs from Clark-2 Model 44 

Table 4.13 Recession Coefficient of Tank Model 46 

Table 4.14 Calculation of Flow Hydrograph by Tank model...........'... 49 

Table 4.17 Auto-correlation coefficients ., 50 

Table 5.1 Performance Evaluation Criteria from Various Models 

(Combined Input) 73 

' Table 5.2 Performance Evaluation Criteria from Various Models: 

Category I ..74 

Table 5.3 Performance Evaluation Criteria from Various Models: 

Category II 75 

Table 5.4 Performance Evaluation Criteria from Various Models: 

Category III 76 

Table 5.5 Performance Evaluation Criteria from Various Models: 

Category IV 77 


X 



Chapter 1 


Introduction 

1.1 General 

Rapid economic and population growth, social changes, urbanization and 
industrialization during the last few decades have resulted in increased and diversified 
demand for water. Due to all these effects, water scarcity has been observed in many 
regions of the world, especially developing countries with arid and a semi arid 
climate. On the other extreme, flood losses have been soaring worldwide. Thus at the 
dawn of the 21*' century, we find ourselves facing a formidable challenge of 
managing our water resources in a sustainable manner. Researchers and experts, the 
world over, engaged planning, development and management of water resources 
systems, have used many modeling tools and techniques such as Systems Analysis 
and Simulation, Geographical Information Systems (GIS), Artificial Neural Networks 
(ANN), Conceptual models. Fuzzy Logic, Decision Support Systems (DSS) and 
Expert Systems for modeling and management of water resources system. The 
continuing evolution of computer science and information technology has helped in 
developing more sophisticated modeling tools for tackling problems of water 
resources. 

1.2 Statement of Problem 

The focus of present research work is on the modeling of the rainfall-runoff process 
using a variety of conventional and new techniques. The rainfall-runoff (R-R) process 
is a highly nonlinear, time varying, spatially distributed, complex and dynamic 
process which is not easily understood and is difficult to model [Singh, 1964; 
Kulandaiswamy and Subramanian, 1967; Chiu and Huang, 1970; Pilgrim, 1976\. 
Two major approaches for modeling the R-R process are normally employed: 
"conceptual" (physical) modeling and "systems-theoretic" modeling (sometimes 
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referred to as "black box") [Amorocho and Hart, 1964; O’Connell and Clarke, 1981; 
Sorooshian, 1983; Gupta, 1984; Young and Wallis, 1985; Singh, 1988; Duan, 1991]. 

R-R modeling attempts to establish a link between rainfall measured over the 
catchment and runoff (streamflow) observed at the outlet of the catchment. The R-R 
models can be implemented in an operational context for many applications of 
hydrological engineering (e.g. water works design) and water resources management 
(e.g. floods or low flows forecasting, reservoir management). R-R models are tools 
increasingly commonly used in flood forecasting applications in an operational 
context and links may be made with the analyses carried out on that topic (e.g. 
Bumash, 1995; Refsgaard, 1997; Tingsanchali and Gautam, 2000; Yang and Michel, 
2000 ). 

R-R modeling can be divided into two parts, continuous modeling and event-based 
modeling. Continuous modeling is generally used for water supply management 
system; irrigation purposes, power generation, recreation, and fish and wild life 
propagation. Whereas event-based modeling is generally useful for flood forecasting, 
construction of hydraulic structure etc. This thesis is focused on event-based. R-R 
modeling using both conceptual and systems-theoretic approach. 

The conceptual rainfall-runoff (CRR) models are designed to approximate within their 
structures (in some physically realistic manner) the general internal sub-processes and 
physical mechanisms, which govern the hydrologic cycle. CRR models usually 
incorporate simplified forms of physical laws and are generally nonlinear, time- 
invariant, and deterministic, with parameters that are representative of watershed 
characteristics. Until recently, for practical reasons (data availability, calibration 
problems, etc.) most conceptual watershed models assumed lumped representations of 
the parameters. Conceptual watershed models are generally reported to be reliable in 
forecasting the most important features of the hydrograph, such as the beginning of 
the rising limb, the time and the height of the peak, and volume of flow. However, the 
implementation and calibration of such a model can typically present various 
difficulties requiring sophisticated mathematical tools, significant amounts of 
calibration data, and some degree of expertise and experience with the problem 
domain. While conceptual models are of importance in the understanding of 
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hydrologic processes, there are many practical situations, such as streamflow 
forecasting, where the main concern is with making accurate predictions at specific 
watershed locations. In such a situation, a systems-theoretic model may be preferred 
over a conceptual model. In the systems-theoretic approach, mathematical models are 
used to identify a. direct mapping between the inputs and outputs without detailed 
consideration of the internal structure of the physical processes. In the class of 
systems- theoretic models, a method to predict the runoff response of the watershed 
on the basis of known meteorological and hydrologic time series could be based on 
the application of ANN. The success with which ANNs have been used to model 
dynamic systems in science and engineering suggests that the ANN approach may 
prove to be an effective and efficient way to model the R-R process in situations 
where explicit knowledge of the internal hydrologic sub-processes is not required. 
Due to this, ANNs have been proposed as efficient tools for modeling and prediction 
in hydrology, as black- box models (Bishop 1994). 


There are numerous applications of ANNs for streamflow forecasting. Kang et al. 
[1993] developed ANNs for daily and hourly streamflow forecasting by selecting 
from among four pre-specified network structures. Kuo-lin Hsu et al. [1993] 
demonstrated the nonlinear ANN model approach to provide better representation of 
the rainfall-runoff relationship of the medium size Leaf River basin near Collins, 
Mississippi. Halff et al. (1993) designed a three-layer feed forward ANN using the 
observed rainfoll hyetographs as inputs and hydrographs recorded by the U.S. 
Geological Survey (USGS) at Bellvue, Washington, as outputs. Some other examples 
of ANN in rainfall-runoff modeling include, Minns and Hall (1996); Shamseldin 
(1997); Dawson and Wilby (1998); Compo et .al. (1999); and Zhang and Govindaraju 
(2000). Most of the ANN applications above concentrate on continuous modeling of 
the R-R process. Jain and Indurthy (2003) modeled event-based R-R process in the 
Salado Creek, San Antonio, Texas using ANNs. However, that study was conducted 
on a small watershed using rainfall and runoff data at 5-minute time interval. 

Although many ANN applications to hydrologic engineering have been reported, 
there are certain aspects in ANN model development for which no certain guidelines 
are available. Most ANN application in engineering employ back-propagation training 
algorithm, a trial and error method of considering different number of hidden neurons. 
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constant value of acceptable error for training purposes, constant value of learning 
rate and momentum correction factors. The back-propagation training algorithm has 
been reported to be inefficient in capturing complex R-R relationship and it does not 
guarantee the global optimum solution. The ANN models are normally developed by 
considering a limited number of training eiTor statistics, which may not be able to 
ensure against under- and/or over-training of the network. Some of the aspects will be 
exposed in this study to develop definite guidelines while training of an ANN model 
for a particular data set. 

1.3 Research Objectives 

The objectives of present research effort are: 

1 . Evaluate the efficiency of conceptual techniques to model the event-based 
R-R process for large watersheds. 

2. Evaluate the efficiency of ANN techniques in developing event-based R-R 
models for large watersheds. 

3. Attempt to develop guidelines and methodologies for the proper training of 
the ANN models. 

4. Compare the performance of various models developed using a wide 
variety of standard statistical performance evaluation criteria. 

Many conceptual models and various architecture of ANN model will be explored. 
The rainfall and runoff data from the Kentucky River basin will be employed to 
develop all the models. A wide variety of standard statistical performance evaluation 
parameters will be employed to evaluate the performance of the developed models. 

1.4 Organization of the Thesis 

This thesis contains six chapters including this chapter. First chapter deals with brief 
description of the CRR models and ANN models for rainfall-runoff modeling 
techniques, statement of problem, objective and organization of thesis. Chapter 2 
provides a brief literature review on both CRR based and ANN based R-R models. 
Chapter 3 presents the description of ANN. Chapter 4 describes the methodology of 
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models developed in this study. Chapter 5 includes the results and discussions. 
Finally, chapter 6 presents summary, conclusion and recommendations for future 
study in this area. References and appendices are provided at the end. 



Chapter 2 


Literature Review 

2.1 General 

A lot of work has been reported on rainfall-runoff modeling. Based on the problem 
focused in this study, the literature review has been divided into two parts: CRR 
models and ANN models. 

2.2 CRR Models 

CRR models are designed to approximate within their structures (in some physically 
realistic manner) the general internal sub-processes and physical mechanisms, which 
govern the hydrologic cycle. CRR models usually incorporate simplified forms of 
physical laws and are generally nonlinear, time-invariant, and deterministic, with 
parameters that are representative of watershed characteristics. Until recently, for 
practical reasons (data availability, calibration problems, etc.) conceptual watershed 
models assumed lumped representations of the parameters. 

The unit hydrograph (UH), a method for estimating storm runoff, was first proposed 
by Sherman in 1932 (Chow et al. 1988), and since then has been used as a key 
concept. The UH was defined as the watershed response to a unit depth of excess 
rainfall, uniformly distributed over the entire watershed and applied at a constant rate 
for a given period. In 1938, after studying watersheds in the Appalachian mountains 
of the United States, Snyder proposed relations between some of the characteristics of 
the unit hydrograph, i.e. peak flow, lag time, base time, and width (in units of time) at 
50% and 75% of the peak flow. Snyder’s method was enhanced with the 
regionalization of the watershed parameters developed in 1977 by Espey, Altman and 
Graves. A significant contribution to the unit hydrograph theory was given by Clark 
(1945), who proposed a unit hydrograph which was the result of a combination of a 
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pure translation routing process followed by a pure storage routing process. Although 
Clark did not develop a spatially distributed analysis, the translation part of the 
routing was based on the time-area diagram of the watershed. The storage part 
consisted of routing the response of the translation through a single linear reservoir 
located at the watershed outlet. The detention time of the reservoir was selected in 
order to reproduce the falling limb of the observed hydrographs. Note that the actual 
travel time of a water particle, according to this approach, is the travel time given by 
the time-area diagram plus the detention time of the reservoir, which is somewhat 
inconsistent. Some years later, Nash (1957) proposed a unit hydrograph equation, 
which is a gamma distribution, i.e. the response of a cascade of identical linear 
reservoirs to a unit impulse. It is important to notice that the method proposed by 
Nash did not model the watershed itself, and was just a fitting technique based on the 
first and second moments of the calculated and observed hydrographs. 

The 1960s brought the introduction of computers into hydrological modeling enabling 
complex water problems to be simulated as complete systems. The first 
comprehensive hydrologic computer model, the Stanford Watershed Model (SWM), 
was developed at Stanford University (Crawford and Linsley, 1966). SWM was 
capable of simulating R-R process at both continuous and event basis. In the late 
1960s, HEC-l was developed by the Hydrological Engineering Center, U. S. Army 
Corps of Engineers to simulate surface runoff processes from precipitation and 
snowmelt over a watershed. Varying rainfall distributions and total storm rainfall was 
applied to different portions of a watershed. HEC-l simulates surface runoff response 
of a river basin to precipitation by representing the basin as a system of hydrologic 
and hydraulic components. While such models ignore the spatially distributed, time 
varying, and stochastic properties of the R-R process, they attempt to incorporate 
realistic representations of the major nonlinearities inherent in the R-R relationships. 
Real-time forecasting rainfall-runoff models were developed in the late 1970s and 
1980s. In 1972, the Soil Conservation Service (SCS) of the US Department of 
Agriculture (USDA) proposed a unit hydrograph model based on a single parameter: 
the lag time between the center of mass of the excess precipitation hyetograpH and the 
peak of the unit hydrograph. The shape of the hydrograph is given by an average pre- 
computed dimensionless unit hydrograph curve or, as a simplification, by triangular 
dimensionless unit hydrograph (Chow et al. 1988). Among the more widely used and 



reported lumped parameter watershed models are the Sacramento Soil Moisture 
Accounting (SAC-SMA) Model of the U.S. National Weather Service [Bumash et al, 
1973]. 

Yet, studying the storm rainfall-runoff relation involves much more than studying the 
unit hydrograph. Consequently, in trying to relax the unit hydrograph assumptions of 
uniform and constant rainfall, and to account for spatial variability of the catchment, 
considerable research has been done in recent years, and many articles dealing with 
these topics can be found in the literature. Sugawara (1972) proposed catchment 
response be simulated using conceptual non-linear water tank model. He attempted to 
simulate the rainfall-runoff process by combining several water tanks. His model had 
one .or two outlets on the sidewall of the'water tank to simulate runoff discharge. He 
succeeded in predicting runoff discharge from rainfall fairly accurately. Hino and 
Hasebe (1984) separated the runoff discharge into three components (overland flow, 
subsurface flow and groundwater flow), predicted each component, and obtained 
good results from conceptual water tank model. Mizumura and Chiu (1985) applied 
the tank model to predict combined runoff of snowmelt and rainfall. They also 
developed this rainfall-runoff model using recession curves. 

Mizumura (1995) proposed a simple water tank model to simulate the rainfall-runoff 
process. The model was applied to two catchments, one in Japan and other in US for 
predicting flood response to rainfall events. The results found from the tank model 
were in good agreement with the observed data. 

Basha (2000) proposed a simple nonlinear conceptual model for simulating R-R 
phenomenon in a catchment. The model is based on the concept of a nonlinear storage 
outflow relationship and assumed a nonlinear loss model. The model was applied to a 
very small semi-arid watershed with considerable success. The proposed model can 
be considered as an improvement over the linear reservoir theory and an attractive and 
useful substitute to the linear reservoir components in more complex conceptual R-R 
models. 
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Conceptual watershed models are generally reported to be reliable in forecasting the 
most important features of the hydrograph, such as the beginning of the rising limb, 
the time and the height of the peak, and volume of flow. However, the 
implementation and calibration of such a model can t 3 rpically present various 
difficulties, requiring sophisticated mathematical tools [Sorooshian et al, 1993], 
significant amounts of calibration data [Yapo et al, 1995], and some degree of 
expertise and experience with the model. Alternatively many researchers have 
focused their attention in the ANN technique to develop hydrologic models. 

2.3 ANN Models 

The use of ANN technique in water resources and streamflow prediction is relatively 
new. ANNs are supposed to possess the capability to reproduce the unknown 
relationship existing between a set of input variables (e.g., rainfall) of the system and 
one or more output variables (e.g., runoff) (Chakraborty et al- 1992). Many studies in 
which ANN models have been applied to problems involving runoff forecasting and 
weather prediction have been reported in the literature. Kang et al [1993] developed 
ANNs for daily and hourly streamflow forecasting by selecting from among four pre- 
specified network structures. The neural network approach was applied to the flow 
prediction of the Hurm River at the Dexter sampling station, Michigan (Karunanidhi 
et al 1994). Empirical comparisons were performed between the predictive capability 
of the neural network models and the most commonly used analytic nonlinear power 
model in terms of accuracy and convenience of use. Preliminary results were found to 
be quite encouraging. Markus et al. (1995) used ANNs with the back propagation 
algorithm to predict monthly river flows at the Del-Norte gauging Station in the Rio 
Grande Basin in south Colorado. The results indicated that ANNs did a good job of 
predicting stream flows. 

Hjelmfelt and Wang (1993a-c) developed a neural network based on the unit 
hydrograph theory. Using linear superposition, a composite runoff hydrograph for a 
watershed was developed by appropriate summation of unit hydrograph ordinates 
and runoff excesses. Rainfall and runoff data from 24 large storm events were chosen 
from the Goodwater Creek watershed (12.2 km^) in central Missouri to train and test 
the ANN. The resulting network was shown to reproduce the unit hydrograph better 
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than the one obtained through the standard gamma function representation. In a later 
study, Hjelmfelt and Wang (1996) compared this method with a regular three-layered 
ANN with back-propagation. They concluded that a regular network could not 
reproduce the unit hydrograph very well and was more susceptible to noise than a 
network whose architecture was more suited for unit hydrograph computations. In an 
application, using two neural networks, Zhu et al. (1994) predicted upper and lower 
bounds on the flood hydrograph in Butter Creek, New York. Smith, and Eli (1995) 
applied a back-propagation neural network model to predict peak discharge and timp. 
to peak over a hypothetical watershed. Data sets for training and validation were 
generated by either a linear or a nonlinear reservoir model. 

Shamseldin (1997) employed neural network technique in the context of R-R 
modeling. The chosen form of network was tested using different types of input 
information. The technique was applied for four different input scenarios in each of 
which some or all of inputs were used. The performance of the technique was 
compared with those of models that utilize similar input information. The result 
suggested that neural network shows considerable promise in the context of R-R 
modeling. 

Dibike and Solomatine (1998) applied ANN approach for downstream flow 
forecasting in the Apure river basin (Venezuela). Two types of ANN architectures, 
namely multi layer-perceptron (MLP) and a radial basis function (RBF) were 
implemented. The performances of these networks were compared with a CRR 
model and they were found to be slightly better for river flow forecasting problem. 
Tokar and Johnson (1999) employed ANN methodology to forecast daily mnoff as a 
function of daily precipitation, temperature, and snowmelt for the Little Patuxent 
River Watershed in Maryland, USA. The sensitivity of the prediction accuracy to the 
content and length of training data was investigated. The ANN rainfall- ranoff model 
compared favorably with results obtained using existing techniques including 
statistical regression and a simple conceptual model. 

ASCE task committee on application of ANNs in hydrology (2000) reported that 
ANNs are robust tools for R-R modeling. After appropriate training, they were able 
to generate satisfactory results for many prediction problems in R-R modeling. 
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Tingsanchali & Gautam (2000) employed a neural network model with a back- 
propagation algorithm for flood forecasting in the tidal reach of the Chao Phraya 
River at Bangkok. The model developed was found to perform very well in both 
calibrations (training) and verification (testing). The efficiency of the forecasting 
model was found to be very satisfactory. Dolling and Varas (2001) employed ANN 
for monthly streamflow forecasting on mountain watersheds. The procedure 
addresses the selection of input variables, the definition of model architecture and the 
strategy of the learning process. 


Jain and Indurthy (2003) investigated the suitability of some deterministic and 
statistical techniques along with the ANN technique to model an event-based rainfall- 
runoff process. Data derived from Salado Creek at Bitters Road, San Antonio, Texas 
were employed. It was found that the ANN models consistently outperformed 
conventional models, barring a few exceptions, and provided a better representation 
of an event-based rainfall-runoff process, in general. However, the duration of the 
rainfall considered was of the order of few minutes only, and the size of the catchment 
was very small. 
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Chapter 3 


Artificial Neural Networks 

3.1 General 

Artificial neural networks (ANNs) are also referred to as neuro-morphic systems, 
parallel distributed processing, neuro-computing, natural intelligent systems and 
machine learning algorithms. An ANN is an information-processing paradigm that is 
inspired by the way biological nervous systems, such as a human brain, process the 
information. The key element of this paradigm is the novel structure of the 
information processing system. ANNs are an attempt at mimicing the patterns of the 
human mind. It is composed of a large number of highly interconnected processing 
elements (neurons) working in unison to solve a specific problem. Thus, it is a 
mathematical model, which has a highly connected structure similar to brain cells. 
They consist of a number of neurons arranged in different layers, an input layer, an 
output layer and one or more hidden layers. The input neurons receive and process the 
input signals and send the output signals to other neurons in the network. The signal 
passing through a neuron is transformed by weights which modify the functions and 
thus the output signal that reaches the following neuron. Modifying the weights for all 
neurons in the network changes the output. Once the architecture of the network is 
defined, weights are calculated to represent the desired output through a learning 
process where the ANN is trained to obtain the expected results. ANNs, like people, 
learn by experience or being presented by manual examples. An ANN is configured 
for a specific application, such as pattern recognition or data classification, through a 
learning process. Advantages of ANNs include adaptive learning, an ability to learn 
how to do tasks based on the data given for training or initial experience, and self- 
organization. 

3.2 An Artificial Neuron 

An artificial neuron is a device with many inputs and one output. The neuron has two 
modes of operation; the training mode and the using mode. In the training mode, the 
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neuron can be trained for particular input patterns. In the using mode, when a taught 
input pattern is detected at the input, its associated output becomes the current output. 
The architecture of an artificial neuron is shown in Figure 3.1. 


neuron in lay^ 



Input Output 

Figure 3.1: An Artificial Neuron 

Note that various inputs to the network are represented by the mathematical symbol, 
Xn. Each of these inputs is multiplied by a connection weight, W^nj- In the simplest 
case, these products are simply summed, fed through a transfer function to generate a 
result, and then output. 

3.3 Architecture of ANNs 

The architecture of an ANN can be of two types i.e. single-layered or multi-layered. 
Based on the training algorithm used to train the network, the architecture can be 
classified as either feed-back network or feed-forward network. Feed-forward ANNs 
allow signals to travel one way only; from input to output. There is no feed-back 
(loops) i.e. the output of any layer does not affect that same layer. Feed-forward 
ANNs tend to be straightforward networks that associate inputs with outputs. They 
are extensively used in pattern recognition. Feed-back networks can have signals 
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traveling in both directions by introducing loops in the network. Feed-back networks 
are very powerful and can get extremely complicated. Feed-back networks are 
dynamic; their state is changing continuously until they reach an equilibrium point. 
They remain at the equilibrium point until the input changes and a new equilibrium 
needs to be found. Feed-back architectures are also referred to as interactive or 
recurrent. 

3.4 Learning 

The learning ability of an ANN is determined by its architecture and by the 
algorithmic chosen for training. The training method usually consists of one of three 
schemes:(l) unsupervised learning: The hidden neurons must find a way to organize 
themselves without help from the outside. In this approach, no sample outputs are 
provided to the network against which it can^ measure its predictive performance for a 
given vector of inputs (2) Reinforcement learning: This method works on 
reinforcement from the outside. The connections among the neurons in the hidden 
layer are randomly arranged, then reshuffled as the network is told how close it is to 
solving the problem. Reinforcement learning is also called supervised learning, 
because it requires a teacher. The teacher may be a training set of data or an observer 
who grades the performance of the network results. Both unsupervised and 
reinforcement learning suffer from relative slowness, inefficiency and reliance on a 
random shuffling to find the proper connection weights. (3) Back propagation: This 
method is proven highly successful in training of multi-layered ANNs. The network is 
not just given reinforcement for how it is doing on a task; information about errors is 
also filtered back through the system and is used to adjust the connections between 
the layers, thus improving performance. 

3.4.1 The Back-Propagation 

Back-Propagation training method is the most widely used method in engineering and 
scientific applications. To begin, the network is initialized; all the connection 
strengths are set randomly within a certain range. The network is then presented with 
some information. Its output is compared with known targets; the responses will 
initially be random. Then working backwards from the output node, each connection 
strength is adjusted so that next time its answer will be closer to the desired one. This 
whole process: input, processing, comparing output with correct answer and adjusting 


14 



connection strengths is called one hack-propagation cycle’, or often just one 
‘iteration’. The network is then presented with another picture and its answer is 
compared with the correct answer, the connection strengths adjusted where needed. 
This process is replaced until all training examples are exhausted. This process of 
training can often take hundreds or thousands of iterations. To test for this, the 
pictures (or whatever input is being used) should be divided into two groups: The 
training set, and the transfer set. The training set is used during back-propagation 
cycles, and the transfer set is used once learning is complete. If the network performs 
as well on the novel transfer stimuli as it did on the training set, then we conclude that 
learning has occurred. 

Backpropagation is a gradient descent algorithm, in which the network weights are 
moved along the negative of the gradient of the performance function. The term 
backpropagation refers to the manner in which the gradient is computed for nonlinear 
multilayer networks. There are a number of variations of the basic algorithm. The 
MATLAB’s Neural Network Toolbox implements a number of these variations. 
Properly trained backpropagation netwoi'ks tend to give reasonable answers when 
presented with inputs that they have never seen. In order to train a neural network to 
perform some task, we must adjust the weights of each unit in such a way that the 
error between the desired output and the actual output is reduced. This process 
requires that the neural network compute the error derivative of the weights (EW). In 
other words, it must calculate how the error changes as each weight is increased or 
decreased slightly. The algorithm computes each EW by first computing the EA, the 
rate at which the error changes as the activity level of a unit is changed. For output 
units, the EA is simply the difference between the actual and the desired output. To 
compute the EA for a hidden unit in the layer just before the output layer, we first 
identify all the weights between that hidden unit and the output units to which it is 
connected. We then multiply those weights by the EAs of those output units and add 
the products. This sum equals the EA for the chosen hidden unit. After calculating all 
the EAs in the hidden layer just before the output layer, we can compute, in similar 
fashion, the EAs for other layers, moving from layer to layer in a direction opposite to 
the way activities propagate through the network. This is what gives back propagation 
its name. Once the EA has been computed for a unit, it is straightforward to compute 
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the EW for each incoming connection of the unit. The EW is the product of the EA 
and the activity through the incoming connection. 


3.5 The ANN Tool Box 


The MATLAB’s ANN toolbox was used in this study to develop the ANN models. A 
brief introduction is provided here. Toolbox is simply a summary of established 
procedures that are known to work well. There are two features of the toolbox that are 
designed to improve network generalization - regularization and early stopping. To 
work in toolbox, first we need to go in MATLAB. A window appears as shown in 
Figure 3.2a. 
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Figure 3.2a: MATLAB window 

Type nntool on MATLAB prompt, a window appears (Figure 3.2b), that is neural 
network tool. It has different space areas such as input, target, output, error, network, 
input delay states, output delay states and have some operating icons. In input space 
area, we can give direct input or we can import data from file, similarly in target. 
Network can also be imported or can be created. New network window will appear as 
shown in Figure 3.2c. There are many types of networks available in toolbox. Feed- 
forward backprop, network is widely used as discussed in previous section. Input 
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range can be provided through input data itself or manually.' This toolbox has a few 
training parameters that need to be specified. 




inputs: 

input2s6 

input2s7 

input2s8 


Networks: 

networki 


Outputs: 




networki ^outputs 


Targets: 
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Errors: 
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Export... 

View 

Delete 




iNttiwurKb uniy 

Initialize... | Sirhulate... j Train... | Adapt... | 



^^11 


Figure 3.2b: Neural Network Tool 

The Levenberg-Marquardt (TRAINLM) function is most widely used training 
function. It was designed to approach second-order training speed without having to 
compute the Hessian matrix. There are two Adaption learning functions, 
LEARNGDM and LEARNGD. LEARNG.DM is Gradient descent with momentum 
weight and bias learning function, it calculates the weight change dW for a given 
neuron from the neuron’s input P and error E, the weight (or bias) W, learning rate 
LR, and momentum constant MC, according to gradient descent with momentum: dW 
= MC*dWprev + (l-MC)*LR*gW. The previous weight change dWprev is stored and 
read from the learning state. LEARNGD is Gradient descent weight and bias learning 
function it calculates the weight change dW for a given neuron from the neuron’s 
input P and error E, and the weight, (or bias) learning rate LR, according to the 
gradient descent: dW = LR*gW. 
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Figure 3.2c: Network type and different function 


Performance. function is MSB that is mean square eiTor: It measures the network’s 
performance according to the mean of squared emors. In addition to it, there are two 
more performance error, MSEREG and SSE. 

'Pransfer Functions: There are three types of transfer functions that are available in the 
MATLAB toolbox. These are LOGSIG, TANSIG, and PURELIN. The function 
LOGSIG generates outputs between 0 and 1 as the neuron’s net input goes from 
negative to positive infinity, Alternatively, multilayer networks may use the tan- 
sigmoid transfer function TANSIG in the applications in which input variables can 
range -oo to + oo. Occasionally, the linear transfer function PURELIN is used in 
backpropagation networks. The three transfer functions are shown in Figure 3,3 a-c, 
respectively. . 

a 


El 

a = lossis(ii) 

Figure 3.3a: Log sigmoid Transfer function 
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a 

a = tanaisf n) 

Figure 3.3b: Tan sigmoid Transfer function 

a 

II 

a = purcliii(ii) 

Figure 3.3c: Linear Transfer function 

After adjusting all the functions when network is viewed, it appears as shown in 
Figure 3.4a. A plot of performance error v/s training cycle can be obtained as shown 
in Figure 3.4c. The trained network can be tested using testing data set. 



Figure 3.4a: 2-5-1 Architecture 
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Figure 3.4b: Training Window 



Figure 3.4c: Plot-Error v/s Training cycle 


Once the ANN architecture has been trained and tested, another one can be tried. 
Finally, the best ANN architecture can be selected from among the various 
architectures investigated. 
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Chapter 4 


Model Development 

4.1 General 

Two types of modeling techniques have been investigated for the purpose of modeling 
the event based R-R process in this study. The first type of models are conceptual 
models and the second type of models are ANN models. Among the conceptual 
category, four different models have been investigated: Unit hydrograph (UH), NASH 
conceptual model, Clark’s conceptual model and a simple non-linear T ank model. In 
ANN models, many different architecture has been investigated. Stream flow and 
rainfall data from Kentucky River Basin, USA were employed to investigate all models 
structures. The performance of all the models was evaluated using six different standard 
statistical parameters. This chapter begins with a brief description of study area and 
data. 

4.2 Study Area and Data 

The data derived fi:om the Kentucky River Basin were employed to calibrate/train and 
validate/test all CRR/ANN modes developed in this study. The Kentucky River basin is 
shown in Figure 4.1. Forty separate counties lie either completely or partially within the 
boundaries of the Kentucky River Basin. The Kentucky River is the sole water supply 
source for several water supply companies of the state. There is a series of fourteen 
Locks and Dams on the Kentucky River, which are owned and operated by the US 
Army Corps of Engineers. The drainage area of the Kentucky River at Lock and Dam 
10 (LD 10) near Winchester is 10244 km^. The data used in this study include average 
daily stream flow (cumec) (US West, 1989a), and daily (average of five stations) total 
rainfall (mm). The five stations are; Lexington airport, Heidelburg, Hyden, Jackson and 
Manchester (Figure 4.1). Out of all the data available data corresponding to 13 isolated 


21 



storms were chosen to develop all models in this study. Out of these 13 storms, data for 
storm 1 (SI) to storm 6 (S6) were used to calibrate/train all CRR/ANN models and the 
data of the remaining storms, storm 7 (S7) to 13 (S13) were used for validation/testing 
of all CRR/ANN models. 



Figure 4. 1: Kentucky River Basin 


LD 10 - Lock and Dam 10 

KRWTP = Kentucky River Water Treatment Plant 

RRWTP = Richmond Road Water Treatment Plant 

These 13 storms were chosen from different years and these were having different 
characteristics i.e. some have low stream flow, some have medium stream flow, and 
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some have high stream flow. Some example storms are graphically presented in Figure 
4.2 through Figure. 4.5. 


Time (days). 



Time (days) 

Fig.urc 4.2: Rainfall and streamflow for the Storm SI (May 1960) 


Time (days) 



Figure 4.3: Rainfall and streamflow for the Storm S2 (May 1961) 
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Time (day) 



TimeTti! 


Figure 4.4: Rainfall and streamflow for the Storm SIO (April 1970) 


Time (day) 



Time (day) 


Figure 4.5: Rainfall and streamflow for the Storm S13 (May 1972) 

The model development was carried out using_two different approaches in this study. 
The first approach used rainfall and runoff data combined from the calibration storms 
(SI to S6). In the second approach the data were divided into different categories based 
on the peak maj^ iude. 
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Effective Rainfall Effective Rainfall 
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Figure 6: Relationship Between Effective Rainfall and Peak Flow 
from Calibration Storms 



Figure 7 : Relationship Between Effective Rainfall and Peak Flow 

from All 13 Storms 
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This was based on the concept that the rainfall-runoff relationships are different for 
different magnitude storms. In order to analyze this, the effective rainfall and flow from 
various storms were plotted. These plots for calibration data set and whole data set are 
shown in Figure 4.6 and 4.7 respectively. 

Thus based on peak flow, all 13 storms were divided into 4 categories. Category-I: 
Storms having peak flow value between 0 to 200 cumec Category-II; Storms having 
peak flow value between 200 to 400 cumec Category-Ill: Storms having peak flow 
value between 400 to 600 cumec Category-IV: Storms having peak flow value between 
600 to 1050 cumec. Each category consisting of at least 1 calibration/training storm and 
1 validation/testing storm. For model development following procedure has been 
applied in this study. 

4.3 Model development 

The development of a mathematical model involves representing any physical system 
using mathematical equations that can be solved either analytically or numerically. 
There can be many types of models that can be tried before selecting the final one. 
Once a type of model is chosen, the numerical data can be used to determine the 
parameters of the model. This process is called calibration of the model, in ANN 
terminology it is also known as training. Once a model has been calibrated/trained, it 
needs to be - validated/tested on the data it has not seen before. Then certain standard 
statistical performance evaluation criteria can be used to evaluate the performance of 
the model. This process of calibration, validation and evaluation of model performance 
can be carried out for different types of models. A comparison of various types of 
models can then be made in terms of statistical performance indices to select the best 
model. Finally, the selected model can be employed for forecasting and embedding in 
the general water-resources application. This process of model development was 
followed in this study and is depicted in Figure 4.8. 
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Figure 4.8: Model development procedure 


27 

















Two types of models were developed in this study: conceptual models and ANN 
models. Among the conceptual model category four different models are being 
investigated: Unit hydrograph, NASH’s conceptual model, Clark’s conceptual model 
and a non-linear Tank model. Various ANN architectures were investigated before 
selecting the final one. The development of all the model structures investigated in this 
study is presented next. 


4.3.1 Unit Hydrograph Model 

A hydrograph is a plot of discharge v/s. time at any point of interest in a watershed, 
usually its outlet. Hydrographs are the measures of a watershed's response to rainfall 
events. The initial precipitation does not contribute directly to flow at the outlet, and 
instead it is stored or absorbed. This is termed the initial abstraction. Direct runoff is 
that portion of the precipitation that moves directly into the channel. Although several 
types of rainfall-runoff models are available, unit hydrograph (UH) model still remains 
a useful tool for flood estimation in many basins where recorded hydrological data are 
not sufficient to support distributed rainfall-runoff models. The UH is a powerful 
conceptual tool for describing rainfall-runoff relationships. This method was first 
suggested by Sherman in 1932 and has undergone many refinements since then. 

A UH is defined as the hydrograph of direct runoff resulting from one unit depth (1 cm) 
of rainfall excess occurring uniformly over the basin at a uniform rate for a specified 
duration (D-hours). In using the deterministic UH model, one can first develop a UH 
using the rainfall and runoff data of an isolated rainfall event; and then use the 
developed UH to calculate the runoff hydrograph for another given storm having same 
duration. The UH of a particular duration can be converted to another UH of a different 
duration using either superimposition or the S-curve method, if needed. DRH is direct 
runoff hydrograph obtained after the base flow separation. The haseflow is delayed 
flow that reaches a stream essentially as groundwater flow. 
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A UH has three functions (Figure 4.9), input function, transfer function, and output 
function. Input function is nothing but the excess rainfall hyetograph. Transfer function 
is the UH itself and output function is the runoff hydrograph. As we apply excess 
rainfall on the UH, we get can direct runoff hydrograph and then depending upon 
topography by adding base flow, we get flow hydrograph. 



Figure 4.9: Rainfall-Runoff Process: UH Conceptual Model 


In the present study the rainfall and runoff data from the 13 storms were used to 
investigate the UH method. The data from six storms (S1-S6) were used to develop a 
one-day UH. Once the six UHs were developed; they were used to predict the 
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validation storm. The infiltration losses were computed by using the index method. 
The calculation of $ - index for any storm is explained next. 

4.3. 1.1 Calculation of $ - Index 

The -index is the constant rate of abstractions (mm/day) that will yield an excess 
rainfall hyetograph (ERH) with total depth equal to direct runoff over the watershed. 

ERH can be obtained by subtracting initial loss and infiltration losses from rainfall 
hyetograph (plot of rainfall intensity v/s time). The value of <zi -index is determined by 
picking a time interval At , judging the number of interval M of rainfall that actually 
contribute to direct runoff, subtracting (p At from the observed rainfall from each 
intervals, and adjusting the value of (p and M as necessary so that the direct runoff and 
excess rainfall are equal: 

Where, /?„, is observed rainfall (mm) in time interval At . The values of <p -index for 
the calibration storms (SI to S6) are presented in Table 4.1. 


Table 4.1: Runoff depth and (p index for the calibration storms 


Storm 

SI 

S2 

S3 

S4 

S5 

S6 

Runoff-Depth (mm) 

5.17 

13.87 

1.80 

22.43 

10.22 

1.11 

(p -Index (mm/day) 

25.63 

6.93 

42.4 

19.57 

13.18 

16.9 


In order to determine the (p -index for the validation storms, many regression models 
were investigated. The regression models attempted to relate (p -index with explanatory 
variables such as base flow at the beginning of the storm (Qg), antecedent cumulative 
rainfall for last 2 days (P2), 3 days (P3), 4 days (P4) and 5 days (P5). The plots of <p- 
index and these variables are shown in Figure 4.10 to Figure 4.14. It is clear that (p- 
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index shows dependence on Qg and does not seem to be related to any of the 
cumulative rainfalls. A regression model of the following form was finally developed. 



for Qg < 100 m^/s 

(4. 2) 

</) - iO mm/day 

for Qg> 100 m^/s 

(4.3) 


The value of /5o and ^6, were found to be 35.4759 and -0.2555 respectively. 
Regression equation and its coefficient were estimated using the regression tool in 
EXCEL. 



Figure 4.10: Behaviour of -Index With Ground Water Flow Magnitude 
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Table 4.2: Computation of 24-hour UH by UH theory 



Observed 


Computed 


Time 

Average 

Rainfall 

Streamflow 

ERH 

Baseflow 

DRH 

24-hour 

(days) 

(mm) 

(cumec) 

(mm) 

(cumec) 

(cumec) 

UH 

I 

0 

51.84 

0 

51.84 

0 

0 

2 

0 

46.98 

0 

46.98 

0 

0 

3 

. 9.2 

65.34 

0 

42.58 

22.76 

4.40 

4 

30.8 

167.4 

5.17 

38.59 

128.81 

24.91 

5 

0.6 

233.28 

0 

34.97 

198.31 

38.36 

6 

0 

169.56 

0 

42.02 

127.54 

24.67 

7 

0.4 

125,55 

0 

49.07 

76.48 

14.79 

8 

0.8 

96.66 

0 

56.12 

40.54 

7.84 

9 

0 

81.54 

0 

63.17 

18.37 

3.55 

10 

0 

70.2 

0 

70.22 

0 

0 

11 

0 

61.83 

0 

61.83 

0 ■ 

0 


Table 4.3: Calibrated UHs from UH Model 


Time 

(days) 

UHl 

UH2 

UH3 

UH4 

UHS 

UH6 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

4.4 

10.39 

6.09 

2.88 

4.46 

3.52 

2 

24.91 

25.51 

62.19 

28.32 

29.86 

14.95 

3 

38.36 

31.48 

30.57 

37.74 

36.47 

55.69 ' 

4 

24.67 

22.32 

11.12 

32.92 

23.06 

30.02 

5 

14.79 

14.33 

5.1 

13.4 

13.74 

10.19 

6 

7.84 

9.46 

2.59 

2.82 

7.64. 

3.77 

7 

3.55 

5.05 

1.03 

0.49 

3.29 

0.75 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


UHs obtained using calibration storms are presented in Table 4.3 and Figure 4.15. 
Base flow was separated according to the Fixed Base method while developing the 
UHs. In the fixed base method, the surface runoff is assumed to end at fixed time or 
N days after the hydrograph peak. The base flow before the surface runoff began is 
projected ahead to the time of peak. A straight line is used to connect this projection 
at the peak to the point of recession limb at time N after the peak. 


N = 0.83A°-^ 


(4.4) 
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Error statistics were computed which are discussed in next chapter 

These UHs were validated using the data of the validation storms. Values of runoff 
depth and (?) -indices for validation storms are shown in Table 4.4. 


Table 4.4: Runoff Depth And Index For Validation Storms 


Storm 

S7 

S8 

S9 

SIO 

Sll 

S12 

S13 

Runoff-Depth (mm) 

10.21 

12.33 

12.89 

22.87 

6.66 

7.90 

2.40 

(p -Index (mm/day) 

13.59 

16.27 

14.71 

19.13 

24.94 

13.7 

31.4 


The validation of the developed UHs was done using two different methods. In the 
first method, an average UH was computed using all six calibrated UHs. This can be 
accomplished by simply calculating the average of UH ordinates from UHl to UH6 
for each time step to get an average or combined UH. The ERs from validation 
storms were applied on average UH to compute DRH corresponding to validation 
storms and statistical parameters were calculated. In the second method, the 13 
storms were divided into four different categories based on peak flow as mentioned 
earlier. Then the validation of the UH in each category was done using data from 
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the corresponding category in a similar fashion. The results in terms of various 
performance evaluation measures are presented and discussed in the next chapter. 

4.3.2 NASH’S Conceptual Model 

The NASH conceptual model assumes the catchment to be made up of a series of n 
identical reservoirs each having the same storage constant (K). K is nothing but 
storage time constant and has the dimensions of time. It is approximately equal to 
the time of travel of flood wave through the channel reach. The first reservoir 
receives a unit volume equal to 1 cm of effective rain from the catchment 
instantaneously. This inflow is routed through the first reservoir to get the outflow 
hydrograph. The outflow from the first reservoir is considered as the input to the 
second; the outflow from the second reservoir is the input to the third and so on for 
all n reservoirs. The outflow hydrograph from the n* reservoir was taken as the 
instantaneous unit hydrograph (lUH). An lUH is the limiting case of a UH of zero 
duration. Thus, it is a fictitious, conceptual unit hydrograph, which represents the 
surface runoff from the catchment due to an instantaneous application of the rainfall 
excess volume of 1 cm. 

From the equation of continuity 


i-Q 


dt 


For a linear reservoir S = AQ; and hence 

dt dt 


Substituting in equation (4.5) and rearranging 


dt 


+ Q = l 


(4.5) 


(4.6) 


(4.7) 
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And the solution of differential equation, where Q and I are function of time t, is 


0- 



(4.8) 


For the first reservoir, the input is applied instantaneously. Hence fov t > 0, I = 0. 
Also cit ; = 0, J k!i = instantaneous volume inflow = i cm of effective rain. Hence 
for the first reservoir equation (4.8) becomes, 


Q 


1 



(4.9) 


For the second reservoir 


Q2 = 



(4.10) 


Here, I = input = Q), thus 



(4.11) 


Similarly, the hydrograph of outflow from the n''‘ reservoir Q„ is obtained as 


a, 


1 


1 


(7j-i)!i5:" 


(•''-FfA' 


(4.12) 


As the outflow from the n'*’ reservoir is caused by 1 cm on excess falling 
instantaneously over the catchment equation (4.12) describe the lUH of the 
catchment. Using the notation u(r) to represent the ordinate of the lUH, equation to 
represent the lUH of catchment is written as 
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Figure 4.16: NASH Model- Cascade of « Linear Reservoirs 


m(0 = 


1 1 


r-'e 


■/k 


Equation (4.13) can be rewritten as 


(4.13) 


uit) = 


1 

/sT(n) 


it/KY:' 



(4.14) 


Here, if t is in hours, u(/') will have the dimension of cm/h; K and n are constants 
for the catchment for to be determined by effective rainfall and flood hydrograph 
characteristics of the catchment. In this study for determination of n and 2$: the 
moment method was applied (Subraraanya, 1997). The sample calculation of n and 
K are shown in Table 4.5 and 4.6; and the computed value of n & K are shown in 
Table 4.7. These values of parameters n and K of NASH model, lUH were obtained 
from equation (4.14) and finally, UHs were obtained which are shown in Figure 
4.17 
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Table 4.5: Calculation of Mj^ and M,^ of Nash’s Model 


Time 

(day) 


Interval 

At 

(day) 

Incremental 

Area 

Moment 

Arm 

First 

Moment 

Second 

Moment 

Part (a) 

Second 

Moment 

Part (b) 

Mn 

Mu 

0 

0 

0 

0 

0 

0 

0 

0 

0.5 

0.333 

1 

5.17 

1 

5.17 

0.5 

2.585 

1.29 

0.431 



2 

0 

1 

0 

1.5 

0 

0 

0 



3 

0 

1 

0 

2.5 

0 

0 

0 



4 

0 

1 

0 

3.5 

0 

0 

0 



5 

0 

1 

0 

4.5 

0 

0 

0 



6 

0 

1 

0 

5.5 

0 

0 

0 



7 

0 

1 

0 

6.5 

0 

0 

0 



8 

0 

1 

0 

7.5 

0 

0 

0 



9 

0 

1 

0 

8.5 

0 

0 

0 



10 

0 

1 

0 

9.5 

0 

0 

0 



11 

0 

1 

0 

10.5 

0 

0 

0 





SUM= 

5.17 


2.585 

1.2925 

0.43083 




Table 4.6: Calculation of Mq^ and Mq^ of Nash’s Model 


Time 

(day) 

DRH 

(cumec) 

Average DR 
Rate in At 

(cumec) 


Incremental 

Area 

Moment 

Arm 

First 

Moment 

Second 

Moment 

Part (a) 

Second 

Moment 

Part (b) 



0 

0 

0 

0 

0 

0 

0 

0 

0 



1 

22.76 

11.38 

1.00 

11.38 

0.50 

5.69 

2.84 

0.95 

3.49 

14.42 

2 

128.81 

75.78 

1.00 

75.78 

1.50 

113.67 

170.56 

6,32 



3 

198.31 

163.56 

1.00 

163.56 

2.50 

408.90 

1,022.25 

13.63 



4 

127.54 

162.92 

1.00 

162.92 

3.50 

570.23 

1,995.83 

13.57 



5 

76.48 

102.01 

1.00 

102.01 

4.50 

459.04 

2,065.70 

8.50 1 



6 

40.54 

58.51 

1.00 

58.51 

5.50 

321.80 

1,769.93 

4.87 



7 

18.37 

29.45 

1.00 

29.45 

6.50 

191.45 

1,244.47 

2.45 



8 

0.00 

9.18 

1.00 

9.18 

7.50 

68.88 

516.65 

0.76 



9 

0.00 

0.00 

1.00 

0.00 

8.50 

0.00 

0.00 

0.00 



.10 

0.00 

0.00 

1.00 

0.00 

9.50 

0.00 

0.00 

0.00 



11 

0.00 

0.00 

1.00 

0.00 

10.50 

0.00 

0.00 

0.00 






SUM= 

612.81 


2139.7 

8788.2 

51.06 


|■||| 


1 Mqi- Ma=nK nK= 2.99162 

2 Mq2- Mi 2 =nK*nK +nK*K +2*nK*Mn /<:= 0.71848 

3 n= nK/K n= 4. 16383 
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Table 4.7: Nash’s Model Parameters- n and K 


STORM 

n 

Kiday) 

(p -Index 

(mm/ day) 

Runoff-Depth 

(mm) 

Time base 

(days) 

SI 

4.1638 

0.7185 

25.63 

5.17 

9 

S2 

3.2390 

0.9106 

6.93 

13.87 

9 

S3 

3.1772 

0.6777 

42.4 

1.80 

9 

S4 

5.2290 

0.5363 

19.57 

22.43 

9 

S5 

3.8818 

0.7480 

13.18 

10.22 

9 

S6 

6.0692 

0.4711 

16.9 

1.11 

9 

Average 

4.2933 

0.6770 

20.77 

9.10 

9 



Figure 4.17: UHs for Nash model 


The 24-hour UHs were developed using these lUHs. This procedure is explained in 
a tabular form in Table 4.9. The calculations of UHs from storms S1-S6 are 
presented in Table 4.8 and are shown in Figure 4.17. Once the UHs were developed 
using calibration storms, they were validated using the data from validation storms 
by computing various performance statistics, using two different methods (i.e. 
combines UH and category wise), as explained earlier. The result in terms of 
statistical parameters are presented and discussed in the next chapter. 
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Table 4.8: Calibrated UHs from NASH Model 


Time 

(days) 

UHl 

UH2 

UH3 

UH4 

UHS 

UH6 


0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


1 

0.01 

0.05 

0.03 

0.00 

0.01 

0.00 


2 

1.04 

2.31 

3.32 

0.52 

1.35 

0.26 


3 

8.43 

11.89 

18.63 

7.63 

9.66 

5.99 


4 

22.72 

24.97 

35.36 

25.21 

24.07 

23.78 


5 

31.30 

29.59 

33.02 

36.32 

31.34 

37.63 


6 

26.82 

23.67 

18.42 

28.44 

25.81 

30.39 


7 

16.28 

14.34 

7.09 

14.04 

15.27 

14.56 


8 

7.65 

7.08 

2.08 

4.89 

7.07 

4.67 


9 

2.95 

3.00 

0.50 

1.30 

2.72 

1.10 


10 

0.98 

1.12 

0.10 

0.28 

0.90 

0.20 


11 

0.29 

0.38 

0.02 

0.05 

0.27 

0.03 


Table 4.9: 

Calculation 24-hour UH from NASH Model 


i 

fe=0.431 

n=4.1638 

r (n)=7.39839 Area= 10244 km^ 







u(t) 

u(t) 

24-h 

Time 

t/k 

(t/k), n-1 

(e, -t/k) 

u(t) 

(cumec) 

lagged 

UH 

(day) 




(cm/day) 


by24h 

(m,3/s) 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1.39 

2.85 

0.25 

0.13 

157.84 

0.00 

78.92 

2 

2.78 

25.51 

0.06 

0.30 

351.68 

157.84 

254.76 

3 

4.18 

92.01 

0.02 

0.27 

,315.36 

351.68 

333.52 

4 

5.57 

228.61 

0.00 

0.16 

194.82 

315.36 

255.09 

5 

6.96 

463.13 

0.00 

0.08 

98.12 

194.82 

146.47 

6 

8.35 

824.56 

0.00 

0.04 

43.43 

98.12 

70.78 

7 

9.74 

1,342.86 

0.00 

0.01 

17.59 

43.43 

30.51 

8 

11.13 

2,048.84 

0.00 

0.01 

6.67 

17.59 

12.13 

9 

12.53 

2,974.03 

0.00 

0.00 

2.41 

6.67 

4.54 

10 

13.92 

4,150.63 

0.00 

0.00 

0.84 

2.41 

1.62 

11 

15.31 

5,611.44 

0.00 

0.00 

0.28 

0.84 

0.56 

12 

16.70 

7,389.76 

0.00 

0.00 

0.09 

0.28 

0.19 
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4.3.3 Clark’s Model 


Based on an extension to time-area (TA) concept, Clark proposed a model for R-R 
transformation in watersheds. This method, also known as Time Area Histogram 
method, aims at developing lUH due to an instantaneous application of unit rainfall 
excess over the catchment. Time here refers to the time of concentration. The time of 
concentration tc is the time required for water from the farthest point of the catchment 
to reach the outlet. It represents the maximum time of translation of the surface runoff 
of the catchment. In gauged catchments, the time of concentration can be taken as the 
time travel between the end of the rainfall excess and the point of inflection of the 
resulting surface runoff. In ungauged catchments the empirical formulae can be used 
to estimate fc.The Clark’s method assumes that the rainfall excess first undergoes pure 
translation and then attenuation. It also assumes that a linear reservoir is 
hypothetically available at the outlet to provide the requisite attenuation. The method 
uses time variable isochrones. A line joining the points having equal times of travel is 
called an isochrone or runoff isochrone. 



Area. 



Figure 4.18(a) and (b): Isochrones and TA histogram for a watershed 


The linear reservoir at the outlet, used for attenuation, can mathematically be 
described by equation S = KQ , where K is the storage time constant. The value of the 
K can be estimated by considering the point of inflection (P,) of the surface runoff 
hydrograph. Point of inflection (P,) is the point at which the inflow into the channel 
has ceased and beyond this point the flow is entirely due to withdrawal from the 
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channel storage. Channel storage is the volume formed by an imaginary plane parallel 
to the channel bottom drawn at the outflow section water surface. Continuity equation 
for linear reservoir is 


I-Q = 


dt 


(4.15) 


Beyond the point of inflection, the flow is entirely due to withdrawal from the channel 
storage. So continuity equation becomes 


dt dt 


(4.16) 


Hence 


K = -Qj{dQldt\ 


(4.17) 


Where suffix i refers to the point of inflection, The inflow rate 1 (cumec) between an 
inter-isochrone area (km^) with time interval (days) is 


AreaxlO^^ 

24x3600xAr^ 


(4.18) 


Knowing K of the linear reservoir, the inflow at various times is routed by the 
Muskingum method (Subramannaya-1997) using the following equation 


Ql ~ (' 0-^2 6 - 1^1 


(4.19) 
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But for linear reservoir, x = 0; (where x is weighting factor) therefore Q = C; and Q 2 
can be written as 


Q2=2CJ,+C,Q, 


(4.20) 


Where 

C, = (0.5At, )/{K + 0.5Ar, ) And C^=(K- 0.5 At ^ )[{k + 0.5At, ) 


In this study, two different Clark models were developed, Clark- 1 model and Clark-2 
model. In Clark- 1 model, the whole catchment was considered as one area. Thus for 
this model, time interval (Ate) for inflow rate between an inter-isochrone area is equal 
to the time of concentration tc. Whereas in Clark-2 model catchment was divided into 
two sub- areas, having Ai=3944 km^ and A2=6300 km^. Clark parameters K, Cj and 
C 2 were calculated for both two models. The values of K and Ci, C 2 for the catchment 
employed in this study are shown in Table 4.10. Using the parameters K, Ci and C 2 of 
Clark- 1 and Clark-2 models, lUHs were obtained. The calculation of lUH and Uh for 
a sample storm are shown in Table 4.1 1 


Table 4.10; Clark’s Model Parameters- Ci, C 2 and K 


Clark- 1 ; A tc = 2 days 

Clark-2; Atc = 

1 day. 


K 

Cl 

C 2 

K 

Cl 

C 2 

Storm 

(days) 



BBS 



SI 

' 2.661 

0.2731 


2.661 

0.1582 

0.6836 

S2 

3.85 

0.2062 


3.85 

0.1149 

0.7701 

S3 

1.058 

0.486 



0.321 

0.358 

S4 

8.897 

0.101 


8.897 ! 

0.0532 

0.8936 

S5 

2.556 

0.2813 


2.556 

0.1636 

0.6727 

S6 

1.944 

0.3397 


1.944 

0.2046 

0.5908 

Average 

3.494 

0.2812 

0.4376 

3.494 

0.1693 

0.6615 
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Table 4.11: Calculation of lUH and UH by Clark’s Method 


Time 

(day) 

Area 

(km^) 

h 

(cumec) 

2Cih 

(cumec) 

C 2 Q 1 

Ordinate 

of lUH 

(cumec) 

lUH 

Lagged by 

24 Hr 

24 Hr 

UH 

(cumec) 

0 

0 

0 

0 

0 

0 

0 

0 

1 

3944 

456.48 

144.25 

0.00 

144.25 

0.00 

72.12 

2 

6300 

729.17 

230.42 

98.67 

329.08 

144.25 

236.67 

3 




225.09 

225.09 

329.08 

277.09 

4 




153.96 

153.96 

225.09 

189.53 

5 




105.31 

105.31 

153.96 

129.64 

6 




72.03 

72.03 

105.31 

88.67 

7 




49.27 

49.27 

72.03 

60.65 

8 




33.70 

33.70 

49.27 

41.49 

9 




23.05 

23.05 

33.70 

28.38 

10 




15.77 

15.77 

23.05 

19.41 

11 




10.78 

10.78 

15.77 

13.28 

12 




7.38 

7.38 

10.78 

9.08 

13 




5.05 

5.05 

7.38 

6.21 


Table 4.12: Calibrated UHs from Clark-2 Model 


Time 

(days) 

UHl 

UH2 

UH3 

UH4 

UHS 

UH6 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

3.61 

2.62 

7.33 

1.21 

3.73 

4.67 

2 

15.44 

11.45 

28.98 

5.45 

15.94 

19.56 

3 

25.69 

19.82 

41.11 

9.97 

26.39 

31.14 

4 

23.33 

19.45 

26.42 

10.85 

23.72 

25.86 

5 

15.96 

14.98 

9.46 

9.69 

15.95 

15.28 

6 

10.92 

11.54 

3.39 

8.66 

10.73 

9.03 

7 

7.47 

8.88 

1.21 

7.74 

7.22 

5.33 

8 

5.11 

6.84 

0.43 

6.92 

4.86 

3.15 

9 

3.49 

5.27 

0.16 

6.18 

3.27 

1.86 

10 

2.39 

4.06 

0.06 

5.52 

2.20 

1.10 

11 

1.63 

3.12 

0.02 

4.93 

1.48 

0.65 
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Figure 4.19: UHs for Clark-2 model 


The calibrated UHs are presented in Table 4.12 and shown in Figure 4.19. Once the 
UHs were obtained, they were validated using the rainfall and runoff data of the 
validation storms using combined approach and category wise, as explained earlier. 
The results in terms of various performance statistics are presented and discussed in 
the next chapter. 

4.3.4 Nonlinear Tank Model 

A simple nonlinear water-tank model (due to Mizumura, 1995) was developed to 
simulate the rainfall-runoff process in the Kentucky River basin. Development of this 
model is provided in detail here. The water tank has a horizontal cross-sectional area, 
which is a function of the water depth and which may be determined from the 
recession curves of flood hydrographs. The total vertical depth of the water tank and 
the infiltration velocity, which governs the time lag between rainfall and mnoff, are 
computed from the rising limbs of the flow hydrograph. 

Total water depth has been divided into three components as per Mizumura Tank 
model, for each component the runoff discharge at time, t + t,is computed as: 
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(4.21) 


Q' - + .Qy^Q 

2 + kciAt 


Q = 


2k^R^t + Q{2-k,M) 
2 + k^At 


;Q,>Q>Q, 


(4.22) 



2k^RAt + Q(2-kAt) 

o 6 

2 + k^At 


;Q,^Q 


(4.23) 


4.3.4. 1 Determination of Parameters 


For prediction of runoff discharge parameters ko, ks, kg and Qs, Qg can be obtained 
from recession curves. The recession curves can be plotted on semi logarithmic paper. 
The slope of each curve can be divided into three parts. The average slope of each 
straight line gives mean recession coefficients ko, ks, kg respectively. The two 
locations where slope of the straight line suddenly changes, determines the discharge 
parameters Qs and Qg. It was found that the slope of the recession curve of the ground 
water flow and subsurface flow components were almost constant and uniform for 
different floods. But the recession curves of the overland flows fluctuated. This shows 
that the overland flow component was produced from the different portions of the 
catchment during different rainfall events. The values of Tank model parameters 
obtained from the six storms are presented in Table 4. 13. 


Table 4.13: Recession Coefficient of non-linear Tank Model 


STORM 

SI 

S2 

S3 

S4 

S5 

S6 

Average 

ko^ 

0.1386 

0.1003 

0.289 

0.0463 

0.1434 

0.1802 

0.1496 

ks= 

0.0957 

0.0896 

0.1719 

0.1621 

0.0887 

0.1054 

0.1189 

kg= 

0.0551 

0.0734 

0.0746 

0.0101 

0.0507 

0.0453 

0.0515 

Qs= 

169.56 

415.8 

59.4 

936.9 

310.5 

56.16 

324.72 

Qs= 

70.2 

182.25 

8.21 

210.6 

137.16 

21.28 

104.95 


46 












Figure 4.20 Determination of Parameters kg, ks, kg, Qg and Qs 


For the computation of mnoff discharge using equations 4.21 to 4.23, the effective 
rainfall values first need to be calculated. Excess rainfall {R) was calculated using 
Green- Ampt method, which is discussed next. 


4.3.4.2 Calculation of Excess rainfall (/?) 

To obtain exact analytical solution of infiltration equations, Green-Ampt proposed a 
more approximate physical theory. The Green-Ampt equation can be represented by 
the following equations: 


Fit)-(pAe\n\ 


(pt^d 


= Kt 


(4.24) 






1 + 


V 


(pA6 ^ 


(4.25) 
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Time (days) 

Figure 4.21: Schematic Plot of Recession Curves on Semi logarithmic Graph 

Green- Ampt parameters were taken from an existing study on the same catchment. 
Using these equation and parameters, infiltration at each time step during the storm 
can be calculated. Once the infiltration has been obtained, effective rainfall at each 
time step can be computed. The values of effective rainfall from all 13 storms are 
presented in Table A-1 1 in Appendix. Knowing the parameters of the Tank model and 
effective rainfall at different time-steps, equations 4.21 to 4.23 were used to compute 
the runoff response of the watershed for different storms. The mnoff was computed 
for both calibration and validation storms using the two different approaches 
(combined and category wise) as explained earlier. A sample calculation of 
computation of flow hydrograph by Tank model is shown in Table 4.14. The results in 
terms of performance statistics are presented and discussed in the next chapter. 
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Table 4.14: Calculation of Flow Hydrograph by Tank model 

kg= 0.05514; ko= 0.13855; ks= 0.09575 per day; Qs= 169.56; Qg= 70.2 cumec 


Time 

day 

Hydrograph 

(cumec) 

ER 

2keR*At 

Q(2-ke*At) 

2+kg*At 

Qf 

2ksR*At 

0(2-ks*At) 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

51.84 

0.00 

0.00 

0.00 

2.06 

0.59 

0.00 

0.00 

2 

46.98 

0.00 

0.00 

1.15 

2.06 

0.56 

0.00 

0.00 

3 

65.34 

0.00 

0.00 

1.09 

2.06 

0.53 

0.00 

0.00 

4 

167.40 

5.71 

0.63 

1.03 

2.06 

0.81 

1.09 

0.00 

5 

233.28 

0.00 

0.00 

1.57 

2.06 

0.76 

0.00 

0.99 

6 

169.56 

0.00 

0.00 

1.49 

2.06 

0.72 

0.00 

0.90 

7 

125.55 

0.00 

0.00 

1.41 

2.06 

0.68 

0.00 

0.82 

8 

96.66 

0.00 

0.00 

1.33 

2.06 

0.65 

0.00 

0.75 

9 

81.54 

0.00 

0.00 

1.26 

2.06 

0.61 

0.00 

0.68 

10 

70.20 

0.00 

0.00 

1.19 

2.06 

0.58 

0.00 

0.62 

11 

61.83 

0.00 

0.00 

1.13 

2.06 

0.55 

0.00 

0.00 


2+ks*At 

Qs 

2koR*At 

Q(2-k„=^At) 

2+ko*At 

Qo 

Qo+Qs+Qg 

(cumec) 

Computed 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

2.10 

0.00 

0.00 

0.00 

2.14 

0.00 

0.59 

70.20 

70.20 

2.10 

0.00 

0.00 

0.00 

2.14 

0.00 

0.56 

66.43 

66.43 

2.10 

0.00 

0.00 

0.00 

2.14 

0.00 

0.53 

62.87 

62.87 

2.10 

0.52 

1.58 

0.00 

2.14 

0.00 

1.33 

157.68 

157.68 

2.10 

0.47 

0.00 

0.00 

2.14 

0.00 

1.24 

146.89 

146.89 

2.10 

0.43 

0.00 

0.00 

2.14 

0.00 

1.15 

136.89 

136.89 

2.10 

0.39 

0.00 

0.00 

2.14 

0.00 

1.08 

127.62 

127.62 

2.10 

0.36 

0.00 

0.00 

2.14 

0.00 

1.00 

119.02 

119.02 

2.10 

0.32 

0.00 

0.00 

2.14 

0.00 

0.94 

111.04 

111.04 

2.10 

0.00 

0.00 

0.00 

2.14 

0.00 

0.58 

68.83 

68.83 . 

2.10 

0.00 

0.00 

0.00 

2.14 

0.00 

0.55 

65.13 

65.13 
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4.3.5 ANN Model 


The main task in developing an ANN model is the selection of the optimal architecture. 
An optimal architecture may be considered as the one yielding the best performance in 
terms of error minimization, while retaining a simple and compact ANN structure. A trial 
and error procedure is generally applied to decide on the optimal architecture. The first 
step in the development of an ANN model is the selection of input and output variables. 

4.3.5.1 Selection of Input and Outputs 

The development of an ANN model involves the selection of the number of input and 
hidden neurons. The identification of input and output neurons can be determined using 
trial and error procedure. In present study, ANN architecture of single hidden layer was 
employed. The output layer consists of a single neuron representing the flow to be 
modelled Q(t) . The neurons were selected based on the auto-correlation analysis of the 
flow data. In a rainfall-runoff modeling, the inputs are normally present and past rainfalls 
and past flow values. The results of auto-correlation are presented in Table 4.15 and 
Figure 4.22. It is clear from these results that the flow value for lag one only is significant 
in modeling Qftj. 


Table 4.15: Auto-correlation coefGcients 


Lag 

(day) 

SI 

S2 

S3 

S4 

S5 

S6 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

0 

1 

1 

1 

1 

1 

1 

1 

0.6243 

0.6525 

0.327 

0.6622 

0.5958 

0.5179 

2 

0.0017 

0.0388 

-0.1445 

0.0599 

-0.0319 

-0.0548 

3 

-0.422 

-0.3945 

-0.3014 

-0,4433 

-0.4987 

-0.3583 

4 

-0.5043 

-0.4884 

-0.1912 

-0.5485 

-0.4331 

-0.4337 

5 

-0.3603 

-0.3958 

-0.1463 

-0.406 

-0.3341 

-0.3017 

6 

-0.2009 

-0.2627 

-0.1328 

-0.1535 

-0.2101 

-0.1561 

7 

0.0212 

-0.0552 

-0.0773 

0.0387 

-0.0103 

-0.0594 

8 

0.1466 

0.1285 

0.0734 

0.1386 

0.1496 

0.1058 
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Lag (day) 

Figure 4.22: Auto-correlation function v/s Lag 


The other input variable in the ANN models were effective rainfalls at two consecutive 
time steps i.e. R( t) and R( t-1). Two different types of ANN models were developed in this 
study. The first ANN model consisting of two input neurons (R(t) and Q(t-1)) while the 
second ANN model consisted of three input neurons (R(t,) R(t-l) and Q(t-1)). Thus the 
two ANN investigated were 2-Nl-l and 3-N2-1. The data need to be encoded and 
normalized before being applied to ANN. In this study, the data were normalized 
between 0 and 1. The values of number of hidden neurons (N1 and N2) were determined 
through training of various architectures and are explained next. 

4.3.5.2 Training of ANN Models 

The purpose of training of an ANN is to determine the set of connection weights and 
thresholds that cause the ANN to estimate outputs that are sufficiently close to target 
values. The data to be employed for training should contain sufficient patterns so that the 
network can mimic the underlying relationship between input and output variables 
adequately. In this study, MATLAB neural network toolbox was used to simulate feed- 
forward back propagation ANN. Following are the MATLAB functions which were used 
in this study: Network type - Feed-forward back propagation. Input ranges: 0-1, Training 

urn Ttn 

^pr^ff5T 
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function - TRAINLM, Adaption learning function - LEARNGDM, Performance 
function - MSB, Transfer function - LOGSIG. 

Most of the ANN applications reported so far employ the trial and error method for 
training the ANN models and determine the optimal architecture for a given data set. 
However, there are many practical difficulties in the trial and error method as there are no 
definite guidelines for the same. For example, there are no definite guidelines regarding 
(a) what should be the acceptable level of global error at the output layer for a particular 
architecture, (b) when the training should be stopped to prevent under-training and over- 
training. In the absence of any definite guidelines, the training process adopted may 
result in (a) an ANN model that is not the best for the given data set, (b) an ANN model 
that is highly sensitive to the initial solution (initial random weights), and (c) an over- 
trained or an under-trained ANN model. An over-trained ANN model would tend to 
perform better during training than in testing, and the trend will be reverse for an under- 
trained network. There are two approaches that are normally used to train the ANN 
models. In the first approach, the total data set is first divided into two parts, training set 
and testing set. The ANN model is trained on the training set and validated on the testing 
data set. However, it is not clear when to stop the training in this approach. Normally, if 
the global error becomes less than or equal to a pre-decided acceptable level of global 
error, then the training is stopped. Another stopping criteria used is the maximum 
number of iterations up to which the program is allowed to train the ANN. In another 
approach, the total data set is divided into three sets, a training set, a validation set, and a 
testing set. The training and validation sets are used to train the ANN model. If we 
consider the global error at the output layer (normally, MSB) as a function of the training 
iterations, it is observed to be a continuously decreasing curve during training. However, 
this curve can start increasing after certain number of iterations for the testing and/or 
validation data set. Therefore, the training is stopped at the number of iterations at which 
the global error starts increasing during the validation set. Although this is expected to 
ensure that the error statistics during training and validation sets will be compatible, it 
may not be able to ensure that the error statistics during the testing data set will also be of 
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function - TRAINLM, Adaption learning function - LEARNGDM, Performance 
function - MSE, Transfer function - LOGSIG. 


Most of the ANN applications reported so far employ the trial and error method for 
training the ANN models and determine the optimal architecture for a given data set. 
However, there are many practical difficulties in the trial and error method as there are no 
definite guidelines for the same. For example, there are no definite guidelines regarding 
(a) what should be the acceptable level of global error at the output layer for a particular 
architecture, (b) when the training should be stopped to prevent under-training and over- 
training. In the absence of any definite guidelines, the training process adopted may 
result in (a) an ANN model that is not the best for the given data set, (b) an ANN model 
that is highly sensitive to the initial solution (initial random weights), and (c) an over- 
trained or an under-trained ANN model. An over-trained ANN model would tend to 
perform better during training than in testing, and the trend will be reverse for an under- 
trained network. There are two approaches that are normally used to train the ANN 
models. In the first approach, the total data set is first divided into two parts, training set 
and testing set. The ANN model is trained on the training set and validated on the testing 
data set. However, it is not clear when to stop the training in this approach. Normally, if 
the global error becomes less than or equal to a pre-decided acceptable level of global 
error, then the training is stopped. Another stopping criteria used is the maximum 
number of iterations up to which the program is allowed to train the ANN. In another 
approach, the total data set is divided into three sets, a training set, a validation set, and a 
testing set. The training and validation sets are used to train the ANN model. If we 
consider the global error at the output layer (normally, MSE) as a function of the training 
iterations, it is observed to be a continuously decreasing curve during training. However, 
this curve can start increasing after certain number of iterations for the testing and/or 
validation data set. Therefore, the training is stopped at the number of iterations at which 
the global error starts increasing during the validation set. Although this is expected to 
ensure that the error statistics during training and validation sets will be compatible, it 
may not be able to ensure that the error statistics during the testing data set will also be of 
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acceptable levels of global error. In this study, all the architectures were trained and 
tested using levels of acceptable error of varying 10'\ 10’^ 10■^ 10'^, 10'*, and 10‘‘°. In 
addition, four standard statistical parameters: coefficient of correlation R, Nash-Sutcliffe 
efficiency E, average absolute relative error (AARE), and root mean square error 
(RMSE) (explained in the next chapter), were computed during both training and testing 
data sets for different levels of acceptable global error. Then the difference between the 
training and testing error statistics (AR, AE, AAARE, and ARMSE) were plotted against 
the acceptable levels of global error for each of the architecture investigated. The level of 
acceptable global error corresponding to minimum differences between training and 
testing errors was selected as the best model for a particular architecture. For example, 
the plots of global error v/s differences in various error statistics from 2-5-1 ANN 
architecture are shown in Figure 4.23 and 4.24. It is clear from these figures that for the 
2-5-1 ANN architecture for the rainfall-runoff problem under consideration, an 
acceptable error level of 10“^ gives the least difference between training and testing 
statistics. Thus, this level of acceptable global error at the output layer of a 2-5-1 ANN 
model for the data set given would result in the best performance that would avoid any 
under-training or over-training of the network. It is to be noted that one needs to look at 
the actual error levels also (and not just the Aerror) to select the best acceptable level of 
global error to be employed during training. For example, it is clear firom these graphs 
that for an acceptable error level of 10‘* also, the Aerror for most of the error statistics are 
minimum but the actual levels of errors (R, E, AARE, and RMSE) were quite high. 

Further, the sensitivity of the initial solution on the final solution for each of the 
architectures was tested by developing fifteen different models for each of the 
architectures having different initial solutions. The differences in training and testing 
error statistics from fifteen different ANN models was computed for each of the 
architectures. The ANN architecture with a starting solution resulting in the least 
difference in training and testing error statistics (and actual levels of errors) was then 
selected as the best ANN model for a particular architecture. For example, for the 2-5-1 
ANN architecture, the difference in training and testing error statistics for fifteen 
different ANN models having different initial solutions are shown in Figures 4.25 


54 



through 4.26. The 8‘*’ network was selected as the best 2-5-1 ANN model for the data set 
considered. Similarly, for each architecture, the best network was chosen and that 
network was trained at optimum level of acceptable error for that- particular architecture. 
The results in terms of various standard statistical performance evaluation criteria from 
various architectures are presented in the Appendix-A and those from the best 
architecture are presented and discussed in the next chapter. 

The above procedure was adopted for computing the statistical results during training and 
testing for various architectures with number of hidden neurons ranging from 2 to 1 1 for 
both the ANN models. The architecture with the certain number of hidden neurons 
performing the best in terms of both training and testing results was then selected as the 
best. 



Acceptable level of global error 

Figure 4.23: Plot of Acceptable level of global error and Difference Between 
Training and Testing AARE (training S6) from ANN 2-5-1 Model 
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Delta Error (R & E) 
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Acceptable level of global error 

Figure 4.24: Plot of Acceptable level of global error and Difference Betvpeen 
Training and Testing E (training S6) from ANN 2-5-1 Model 
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Figure 4.25: Difference In Performance During Training And Testing For ANN 2-5- 
1 Model For Different Initial Set Of Weights 
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Delta E 



Delta Error (AARE & RMSE) 



Different Set of Initial Weights' 

Figure 4.26: Difference in Performance during Training and Testing for ANN 2-5-1 
Model For Different Initial Set Of Weights 



Observed Streamflow Qo (cumec) 

Figure 4.27: Scatter Plot of Observed and Difference Between Observed and 
Predicted Streamflow from UH Model During Validation 
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Figure 4.28: Scatter Plot of Observed and Difference Between Observed and 
Predicted Streamflow from NASH Model During Calibration 



Observed Streamflow Qo (cumec) 


Figure 4.29: Scatter Plot of Observed and Difference Between Observed and 
Predicted Streamflow from NASH Model During Validation 
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Observed Streamflow Qo (cumec) 


Figure 4.30: Scatter Plot of Observed and Difference Between Observed and 
Predicted Streamflow from ANN 3-10-1 Model During Training 



Observed Streamflow Qo (cumec) 


Figure 4.31: Scatter Plot of Observed and Difference Between Observed and 
Predicted Streamflow from ANN 3-10-1 Model During Testing 
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Figure 4.32: Behaviour of el With Flow Magnitude from UH Model 



Observed Streamflow (cumec) 

Figure 4.33: Behaviour of el With Flow Magnitude from ANN 2-5-1 Model 
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Figure 4.35: Behaviour of e2 With Flow Magnitude from ANN 2-5-1 Model 
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Mathematically it can be expressed as: 



( 5 . 1 ) 


=■ 


1 ^ / 

^ — I ^ Qo,mean )(lQc,x ~ Q-c,mean ) 


( 5 . 2 ) 



( 5 . 3 ) 


5c = 



( 5 . 4 ) 


Where ^ is correlation between observed and computed, is standard deviation of 

observed, 5^ is standard deviation of computed, is observed streamflow, Q^^ is 

computed streamflow, N is total number of data, Qo,mean is the mean of the observed 
streamflow and Qc,mean is the mean of the computed streamflow. 

(ii) Nash-Sutdiffe Efficiency (E): 

Some time coefficient of correlation may not be enough to quantify the strength of 
relationship. Nash-Sutcliffe efficiency is another statistical parameter normally 
employed. Mathematically it can be represented by the following equation: 




Where 

= E ^Qo,i - Qo.,neanf Initial variancc 


( 5 . 5 ) 
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^2 - '^iQc.i Residual variance 

The evaluation schenie for the Nash/ Sutcliffe efficiency is normally taken as follows: 
Efficiency E goodness of fit 


<0.2 

Insufficient 

0.2 - 0.4 

Satisfactory 

0.4 - 0.6 

Good 

0.6 -0.8 

Very good 

>0.8 

Excellent 


(iii) Average Absolute Relative Error (AARE) 

AARE is the average of the absolute values of the relative errors in forecasting a certain 
number of DRH ordinates. To find AARE, we need to first find the relative error (RE) in 
forecasting a data point. RE is a measure of the error in forecasting a particular variable 
relative to its exact value. Mathematically, following equation can represent it: 


RE{i) = 


Qoj 


(5.6) 


Where, RE(i) is the relative error in forecasting, Qo.i is the observed streamflow and 
Qc,i is the computed streamflow. The relative error RE(i) can be either positive or 
negative. Using absolute relative error, AARE can be evaluated as follows: 


AARE = j'^\RE{i)\ 


(5.7) 


Where, N is the total number of data point computed. It is obvious, lower AARE value 
represents good model performance, and vice versa. 
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(iv) Threshold Statistics (TS) 


It is another important parameter to quantify the performance of a model. It measures the 
model performance at a certain level of absolute relative error (ARE). The threshold 
statistics can be defined as the percentage of data points predicted for which the ARE is 
less than a certain level of relative error (say p%). Mathematically, TSp can be 
represented by, 


TS^ =-^^xl00% 
" N 


(5.8) 


Where np is the number of data points whose ARE is less than p% and N is the total 
number of data points. It is obvious that higher the threshold statistics values the better 
the model performance. In this study, the TS statistics has been calculated at ARE 1%, 
5%, 10%, 25%, 50%, and 100%. 


(v) Normalized Mean Balance Error (NMBE) 

NMBE is the average of measure of under-prediction or over-prediction in estimating the 
runoff discharge. Mathematically, it can be represented by the following equation: 


(a., - a.,) 

NMBE = ^ xl00% 

1 N 


The positive value of NMBE represents over-prediction and vice-versa, 

(vi) Root Mean Square Error (RMSE) 


(5.9) 
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The RMSE is another statistical parameter which is employed to measure the residual 
error from model. Mathematically it is expressed as follows: 


RMSE: 


1 


i=l 


N 


(5.10) 


5.3 Discussion of Results 

The results from different models developed in this study were obtained in two ways. 
First, in terms of six different standard statistical parameters and second in terms of 
scattered plot. Results were calculated for both calibration/training data set and 
validation/testing data set to choose the best model. Few sample results from each type of 
model are presented in this chapter and remaining results are provided in Appendix. The 
results in tenns of the various performance evaluation criteria from all the models 
structures are presented in Appendix-A and those from the best models are presented in 
this chapter. The results in terms of statistical parameters from combined or averaged 
approach are presented in Table 5.1 and those from the four categories are presented in 
Table 5.2 through 5.5, respectively. The discussion of the results has been divided into 
two parts: results for combined input and the results for category- wise input. 

5.3.1 Results for Combined Input 

The results in terms of R, E, AARE, TS, NMBE, and RMSE from various models for 
combined. input during both calibration and validation are presented in Table 5.1. It can 
be observed from Table 5.1 that during calibration/training among the CRR models, the 
UH model performed the best and among the ANN models the 3-10-1 ANN model 
performed the best. The UH model obtained the best values of R, E, AARE, NMBE, and 
RMSE statistics of 0.9850, 0.9701, 11.56%, 0.02, and 36.28, respectively, and the 
performance of the UH model was the best in terms of various TS statistics also. The 
performance of the NASH model was comparable to that of the UH model but the 
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performances of Clark- 1, Clark-2, and the TANK models were only reasonable. Further, 
the performances of the UH and the 3-10-1 ANN model was comparable in terms of R, E, 
NMBE and RMSE statistics but UH model performed better than the 3-10-1 ANN model 
in terms ot TS and AARE statistics. Thus, the UH model can be considered to be the best 
during calibration/training when the rainfall and runoff data combined from all the storms 
were used as input in developing the various models. 


Analyzing the results during validation/testing from Table 5.1, it can be observed that the 
3-10-1 ANN model performed the best in terms of R, E, NMBE, and RMSE statistics 
(0.9610, 0.9233, 0.00, and 55.82, respectively) while the UH model performed the best in 
terms of most of the TS and AARE (17.66%) statistics. However, it is worth mentioning 
that all the flow values predicted from the 3-10-1 ANN model had absolute relative error 
(ARE) of less than 100% (see TSIOO in Table 5.1) while no other model was able to 
achieve this. Also, the 3-10-1 ANN model was completely unbiased towards under- 
prediction or over-prediction (NMBE = 0.0). Thus, it can be said that even through the 
UH model performed marginally better than the 3-10-1 ANN model during 
calibration/training, it was not able to perform better during validation/testing. 
Considering the fact that the performance of the 3-10-1 ANN model was comparable to 
the UH model during calibration/training, it can be said that the 3-10-1 ANN model was 
the best in modeling the event-based rainfall-ranoff process in the Kentucky River basin 
when data from many storms varying in peak magnitudes are combined while developing 
the models. This demonstrates the robustness of the ANN technique in developing the 
complex, dynamic, and highly non-linear rainfall-runoff relationships when presented 
with data of varying magnitudes and nature. 

5.3.2 Results for Category-wise Input 

The rainfall and mnoff data from all the available storms were divided into four 
categories depending upon the peak magnitude of the flow values as discussed earlier. 
All the model structures were investigated using the category-wise data in order to test 
their suitability of various models in modeling the complex rainfall-mnoff process of 
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varying degree and complexity depending upon the peak magnitude. The data from all 
the storms were divided into four categories. The results from the four categories are 
presented in Table 5.2 through Table 5.5, respectively. The discussion of results for each 
category is discussed next. Please note that the calibration results from the UH model are 
not included in these tables as when a UH is applied to the calibration storm it would 
produce the same storm giving the perfect error statistics. 

Results for Category-I 

The results in terms of various performance evaluation criteria from various models 
developed using the data in Category-I are presented in Table 5.2. It is clear from Table 
5.2 that the NASH model performed the best during calibration in terms of all statistical 
parameters (apart from the UH model) except for NMBE, for which the 2-5-1 ANN 
model performed the best. The NASH model obtained the best values of R, E, AARE, 
and RMSE statistics of 0.8779, 0.8707, 7.03%, and 6.62, respectively, during calibration. 
Analyzing the results during validation from Table 5.2, it can be observed that the UH 
model performed the best during validation in terms of most of the performance statistics 
barring a few exceptions. For example, the UH model obtained the best R, AARE, 
NMBE, and RMSE statistics of 0.8871, 12.12%, 0.06, and 16.75, respectively; and the 
NASH model obtained the best E statistic of 0.8451. Although the performances of both 
the ANN models were reasonable, they were not better than the CRR models for 
category-I. Thus, it can be concluded that the performances of the conceptual models 
was better than the ANN models and the UH model can be used to model the event-based 
rainfall-runoff process for Category-I for the data under consideration. 

Results for Category-II 

The results in terms of various performance evaluation criteria from various models 
developed using the data in Category-II are presented in Table 5.3. It can be observed 
from Table 5.3 that during calibration/training, the NASH model performed the best in 
terms of E, AARE, and RMSE statistics (values 0.9650, 7.20%, 10.75, respectively), and 
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most of the TS statistics (apart from the UH model); Clarke-2 model performed the best 
in terms of R (0.9740), and the 2-6-1 ANN model performed the best in terms of NMBE 
statistic (0.87). Thus, there was no clear cut winner among all the models developed for 
Categoi 7 -II during calibration although UH and NASH models may be preferred over the 
other models. Analyzing the results during validation/testing for Category-II from Table 
5.3, it can be observed that the NASH model performed the best in terms of R, E, and 
AARE statistics (values 0.9770, 0.9425, and 8.18%, respectively) and some TS statistics; ■ 
the UH model performed the best in terms of NMBE and RMSE statistics (values -0.11 
and 10.49, respectively), and some of the TS statistics. Again the performance of both 
the ANN models was found to be only moderate in comparison to the CRR models 
during validation also. In summary, it can be said that the conceptual models performed 
better than the ANN models and either UH or NASH model can be used to model the 
event-based rainfall-mnoff process for Category-II. 

Results for Category-Ill 

The results in terms of various performance evaluation criteria from various models 
developed using the data in Category-III are presented in Table 5.4. It can be observed 
from Table 5.4 that during calibration/training, the NASH model performed the best in 
terms of R, E, AARE, NMBE, and RMSE statistics (values 0.9891, 0.9739, 5.89%, -0.12, 
and 18.59, respectively), and all of the TS statistics. Further, the 2-10-1 ANN model was 
the joint best in terms of TSl and TS5 statistics. Clearly, the performance of the NASH 
model was excellent and was the best model during calibration for Category-HI (apart 
from UH). Analyzing the results during validation/testing for Category-HI from Table 
5.4, it can be observed that the UH model consistently outperformed the rest of models in 
terms of all the statistical parameters. It obtained the best values of R, E, AARE, NMBE, 
and RMSE statistics of 0.9902, 0.9790, 3.08%, 0.01, and 15.40, respectively. The next 
best model was 3-3-1 ANN model in terms of AARE and NMBE statistics; and the 
NASH model in terms of R, E, and RMSE statistics. The performance of the ANN 
models was found to be comparable to the other CRR models during both calibration and 
validation for Category-III. In summary, it can be said that the UH model is the most 
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suitable model tor modeling the event-based rainfall-runoff process for Category-Ill; the 
performance of the NASH and 3-3-1 ANN model was the next best; and the performance 
of the other CRR and ANN models was moderate and comparable to each other. 


Results for Category-IV 

The results in terms of various performance evaluation criteria from various models 
developed using the data in Category-IV are presented in Table 5.5. It can be observed 
from Table 5.5 that during calibration/training, the NASH model performed the best in 
terms of R, E, AARE, and NMBE statistics (values 0.9793, 0.9709, 8.09%, and -0.11, 
respectively), and all of the TS statistics; and the 2-3-1 ANN model performed the best in 
terms of the RMSE statistics. Further, the 3-3-1 ANN model was the next best model in 
terms of most of the performance statistics. Clearly, the performance of the UH model 
was excellent and it was the best model during calibration for Category-Ill; NASH model 
was the next best; and the 3-3-1 ANN model was the next best model after UH and 
NASH models. Analyzing the results during validation/testing for Category-III from 
Table 5.5, it can be observed that the UH model consistently outperformed the rest of 
models in terms of all the statistical parameters. It obtained the best values of E, AARE, 
NMBE, and RMSE statistics of 0.8926, 9.17%, 0.00, and 51.52, respectively); and the 
NASH model obtained the best R value of 0.9716. The performance of the ANN models 
was not found to be very good. In summary, it can be said that the UH model is the most 
suitable model for modeling the event-based rainfall-runoff process for Category-IV; and 
the performance of the NASH model was the next best. 

5.3.3 Analysis of Various Performance Statistics 

In most of the CRR and ANN applications, the performance of the models is normally 
evaluated using a few error statistics such as R, MSE, and RMSE, etc. The performance 
of various models developed was analyzed using a wide variety of standard performance 
statistics in this study. The error statistics normally employed, such as R, MSE, and 
RMSE, involve the square of the deviations between predicted and observed values of the 
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variable being modeled. Therefore, such statistics tend to be biased towards the high 
magnitudes of the variable being modeled. As such, these statistics tend to be dominated 
by the high magnitude flows and the low magnitude flows are not properly represented in 
them. The error statistics that are not biased towards either high magnitude or low 
magnitude would be more desirable. The average absolute relative error (AARE) is one 
such statistic that is relative to the observed values and is not biased towards any 
magnitude flow. In order to verify this claim, an analysis of these error statistics is 
needed. Such an analysis was earned out in this study and the findings are reported here. 

It is expected that as the magnitude of flow ( Qo) increases, the difference between the 
observed flow and predicted flow (A{Qo-Qp}) would also increase. This can be easily 
verified by plotting the observed flow against the corresponding difference in the 
observed and predicted flow. The QO versus A{Qo-QpJ plots from UH, NASH, and 
ANN models are shown in Figure 4.27 through 4.31, respectively. These figures clearly 
indicate that as the magnitude of flow increases, the deviation between observed and 
predicted flow would also increase. This would mean that the error statistics such as 
MSE and RMSE would be dominated by the high magnitude flows. However, that would 
not be the case for unbiased statistics such as AARE. In order to verify that the AARE is 
a better performance statistic than RMSE, a further analysis of the results is needed. 
Such an analysis should be able to evaluate the impact of the magnitude of the flow value 
on the error statistic (AARE and RMSE). A performance statistic can be considered a 
better one if it does not vary with the magnitude of the variable being modeled. 

Let us define two additional error statistics el and e2. The error statistic el is nothing but 
the ratio of the absolute relative error in predicting a flow value (ARE(i)) to the overall 
AARE in predicting a total of N flow values. Similarly, the error statistic e2 is defined as 
the ratio of the contribution of an individual data point to the RMSE (RMSE(i)) to the 
overall RMSE in estimating a total of N flow values. Figure 4.32 through Figure 4.33 
show the el versus Qo plot from UH and ANN models; while the Figures 4.34 through 
Figure 4.35 show the e2 versus Qo plots from the UH and ANN models. It can be seen 
from these plots that as the value of the observed flow increases, the contribution to 
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RMSE also increases but the conhibuta to overall AAEE ftom an individual data point 
remains constant. This clearly suggests that AARE is an unbiased performance statistic 
that does not depend upon the magnitude of flow and should be preferred over the 
conventional use of error statistics such as MSE and RMSE etc. 
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Table 5.1; Petformance Evaluation Criteria from Various Models (Combined Input) 

STORM MODEL R E AARE TSl TS5 TSIO TS25 TS50 TSIOO NMBE RMSE 

During Calibration/Training 
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Table 2; Performance Evaluation Criteria from Various Models: Category I 

E AARE TSl TS5 TSIO TS25 TS50 TSIOO NMBE RMSE 
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Table 3: Performance Evaluation Criteria from Various Models: Category II 

E AARE TSl TS5 TSIO TS25 TS50 TSIOO NMBE RMSE 
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Table 4. Performance Evaluation Criteria froni Various Models; Category III 
E AARE TSl TS5 TSIO TS25 TS50 TSIOO 
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Figure 5.1: Scatter Pot for UH-c Model During Calibration 



Figure 5.2: Scatter Pot for UH-c Model During Validation 
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Figure 5.3: Scatter Plot for ANN 2-8-lc During Training 
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Figure 5.4: Scatter Plot for ANN 2-8-lc During Testing 
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Figure 5.5: Scatter Plot for ANN 3-10-lc During Training 



Figure 5.6: Scatter Plot for ANN 3-10-lc During Testing 
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Figure 5.11: Scatter Plot for Category-I storms from ANN 2-5-1 Model During 

Training 



Figure 5.12: Scatter Plot for Category-I storms from ANN 2-5-1 Model During 

Testing 
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Figure 5.13: Scatter Plot for Category-II storms from ANN 2-6-1 Model During 

Training 
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Figure 5.35: Scatter Plot for Category-Ill storms from ANN 2-10-1 Model During 
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Figure 5.16: Scatter Plot for Category-Ill storms from ANN 2-10-1 Model During 
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Figure 5.17: Scatter Plot for Category-IV storms from ANN 3-6-1 Model During 

Training 



Figure 5.18: Scatter Plot for Category-IV storms from ANN 3-6-1 Model During 
' Testing 
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Figure 5.19: During Calibration/Training Combined Input 
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Figure 5.21: Category-I Storm During Calibration/Training 
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Figure 5.22: Category-I Storm During Validation/Testing 
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Figure 5.25: Category-Ill Storm During Calibration/Training 
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Figure 5.26: Figure 2: Category-Ill Storm During Validation/Testing 
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Figure 5.27: Category-IV Storm During Calibration/Training 



Figure 5.28: Figure 2: Category-IV Storm During Validation/Testing 
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Chapter 6 


Summary, Conclusions, and Scope of Future 
Work 

6.1 Summary 


This thesis presents the findings of a study on the evaluation of various techniques 
available for modeling the event-based rainfall-runoff process. The techniques 
'investigated include conceptual techniques and artificial neural network (ANN) 
technique. Among the conceptual techniques, four different conceptual rainfall runoff 
(CRR) model stmctures were investigated: unit hydrograph (UH), NASH’s conceptual 
hydrograph model, Clark’s conceptual hydrograph model, and a non-linear water TANK 
model. Among the ANN models, two types of ANN models were investigated that 
differed in the number of input variables considered to model the complex, dynamic, and 
non-linear event-based rainfall-mnoff process. The rainfall and flow data derived from 
the Kentucky River basin were employed to develop all model structures investigated in 
the present study. The data considered in this study include spatially averaged daily total 
rainfall (mm) taken from five different locations scattered throughout the Kentucky River 
basin, and the daily average flow on Kentucky River near Winchester, Kentucky, USA. 
The total rainfall and average flow data from thirteen isolated storms were available, 
which were divided into two sets: a calibration or training set consisting of data from six 
storms, and a validation or testing set consisting of the data from the remaining storms. 
The whole data set was also sub-divided into four categories depending upon the peak 
magnitude in order to test the hypothesis that the models developed using the 
decoipposed data sets based on physical concepts may perform better than the models 
developed on the whole data set consisting of a wide variation in flow magnitudes and 
hence more complex rainfall-runoff relationships. An attempt was also made to develop 
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guidelines to ensure that the developed ANN models are the best models for a given data 
set while training the various ANN models, and that the developed ANN models are not 
over-trained or under-trained. A wide variety of standard statistical performance 
evaluation measures were employed to evaluate the performances of all the model 
structures investigated. 

6.2 Conclusions 

The results of the evaluation of various techniques for rainfall-runoff modeling carried 
out in this study for the case when the data from many storms are combined into one 
calibration set indicate that the 3-10-1 ANN model was the best in terms of various 
performance statistics. The UH model performed comparable to the best ANN model 
during cahbration/training, the ANN model performed significantly better than the CRR 
models during validation/testing. These results strengthen and verify the earlier findings 
by many researchers that the ANN technique is quite efficient in modeling the complex, 
dynamic, and non-linear rainfall-mnoff process. The massively parallel stmcture of the 
ANN models help them in better generalizing the complex rainfall-runoff relationships 
inherent in a large data set consisting of wide range of flow magnitudes. Since the 
rainfall-runoff relationships in different magnitude events are different, it is better to rely 
on a model structure that can not only provide good efficiency in modeling but also good 
effectiveness in accurately predicting the flow values. 

The results of the evaluation of various techniques for rainfall-runoff modeling carried 
out in this study for the case when the data from many storms are divided into four 
categories depending on the peak magnitude indicate that the conventional modeling 
techniques of UH and NASH are more suitable than the other techniques for modeling 
the event-based rainfall-runoff process in a large watershed. It has been found that (a) the 
performances of the conceptual models were better than the ANN models, and that the 
UH model is the most suitable model for modeling the event-based rainfall-runoff 
process for the data under Category-I; (b) the conceptual models performed better than 
the ANN models for modeling the event-based rainfall-runoff process under Category-II, 
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and the UH and NASH models were found to be better than all other models; (c) the UH 
model was found to be the most suitable model for modeling the event-based rainfall- 
runoff process for Category-Ill; the performance of the NASH and 3-3-1 ANN model 
was the next best, and the performance of the other CRR and ANN models was moderate 
and comparable to each other; and (d) the UH model is the most suitable model for 
modeling the event-based rainfall-runoff process for Category-IV; and the performance 
of the NASH model was the next best in modeling the event-based rainfall-runoff model. 
In summary, it can be said that when the data are divided into different categories based 
on peak flow values, the conceptual methods can be preferred over the ANNs and that the 
UH is still the best method of modeling the event-based rainfall-runoff process. This 
finding is particularly interesting since the UH method is generally considered applicable 
for watersheds of areas up to approximately 5000 km^. However, the drainage area of the 
Kentucky River considered in this study was more than 10,000 km^. 

In this study, an attempt was made to develop new guidelines for training of various 
architectures. It has been found that while training the ANN models, the emphasis should 
be on the evaluation of the models by examining the differences between the training and 
testing error statistics using a wide variety of performance statistics rather than training 
the ANN models using the global error at the output layer during training only. An ANN 
model trained using MSB or similar error statistic during training data set only will tend 
to be either over-trained or under-trained depending upon its architecture. Further, it has 
been found that the performance of the ANN model depends upon the level of acceptable 
global error at the output layer. Instead of using a constant level of acceptable error for 
every architecture during training, it should be made adaptive and an optimum level of 
acceptable global error at the output layer should be identified to prevent over-training 
and/or under-training of the network. Further, the investigation of various error statistics 
considered in this study indicate that the unbiased error statistics such as average absolute 
relative error (AARE) are better indicator of the performance of a model than the error 
statistics normally employed such as MSB and RMSE. This is due to the fact that the 
error statistics such as MSB and RMSE tend to increase with an increase in the 
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magnitude of the variable being modeled, and hence get dominated by the high 
magnitudes of the variables being modeled. 

6.3 Scope for Future Work 

The author would like to mention that no study is complete and there are always scopes 
for improvements and further work. In light of the present research work, the following 
limitations/scope for future works are identified. The conceptual models considered in 
this study are lumped in nature. The performances of the CRR models can be further 
improved by considering the spatial variations in the rainfall, runoff, and other hydrologic 
variables. The infiltration was estimated using the 0-index method in developing the 
UH, NASH, and Clarks’ s models. It would be possible to enhance the performances of 
these models if infiltration is estimated using more accurate methods such as Green- Ampt 
method. It was interesting to find the UH and NASH models performing very well in a 
large watershed such as Kentucky River basin. The suitability of the conceptual 
techniques such as UH and NASH models must be investigated for modeling the event- 
based rainfall-runoff process in large watersheds of varying hydrological and climatic 
conditions. The training of the ANN models was carried out using the back-propagation 
training algorithm, which does not guarantee the global optimum solution and has been 
reported to be inefficient. Although an attempt was made to obtain a solution close to the 
global optimum solution by developing fifteen different ANN models for each 
architecture investigated having different initial weights, the global solution can never be 
guaranteed. It would be interesting to employ other training methods that may improve 
the performance of the ANN models. Finally, the results reported in this study are based 
on a small data set of thirteen storms only; however, the study was limited due to the 
availability of the data. Ideally, it would be better to consider a larger data set and data 
from more watersheds in order to have greater confidence in the findings of such studies. 
It is hoped that the future research efforts will concentrate on some of these directions to 
enhance the performance of the event-based rainfall-runoff models and hence the 
important water resources planning, design, and management activities. 
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Appendix-A 

Performance of Various Models in Terms of Standard 

Statistical Parameters 
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Table A-2; Performance Evaluation Criteria from Various Models (Combined Input) 
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Table A-3; Performance Evaluation Criteria from Various Models: Category I . 

E AARE TSl TS5 TSIO TS25 TS50 TSIOO NMBE RMSE 
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ANN 3-10-1 0.7318 0.2829 35.84 0.0 9.1 18.2 36.4 81.8 100.0 18.74 17.15 

ANN 3-11-1 0.7687 0.7329 102.87 0.0 0.0 9.1 18.2 27.3 45.5 68.54 29.59 


Table A-4: Performance Evaluation Criteria from Various Models; Category I 

E AARE TSl TS5 TSIO TS25 TS50 TSIOO NMBE RMSE 
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Table A-5; Performance Evaluation Criteria from Various Models: Category 11 

model R E AARE TSl TS5 TSIO TS25 TS50 TSIOO NMBE RMSE 

During Caiibration/Training 
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Appendix-B 

Performance of Various Models-Scatter Plots 
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Observed Streamflow (cumec) 


Figure B-1: Scatter Plot for Category-I storms from ANN 3-4-1 Model During 

Training 



Observed Streamflow (cumec) 

Figure B-2: Scatter Plot for Category-I storms from ANN 3-4-1 Model During 

Testing 
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Predicted Streamflow (cuinec) 3 Predicted Stream. 
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Observed Streamflow (cumec) 


B-1: Scatter Plot for Category-I storms from ANN 3-4-1 Model During 

Training 
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Observed Streamflow (cumec) 

Figure B-2: Scatter Plot for Category-I storms from ANN 3-4-1 Model During 

Testing 
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Predicted Streamflow (cumec) 
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FigureB-4: Scatter Plot for Category-II storms from ANN 3-7-1 Model During 

Testing 
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Figure B-5: Scatter Plot for Category-Ill storms from ANN 3-3-1 Model During 

Training 



Figure B-6: Scatter Plot for Category-Ill storms from ANN 3-3-1 Model During 

Testing 
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Figure B-7: Scatter Plot for Category-Ill storms from ANN 3-7-1 Model During 

Training 



Observed Streamflow (cumec) 

Figure B-8: Scatter Plot for Category-in storms from ANN 3-7-1 Model During 

Testing 
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Figure B-9: Scatter Plot for Category-IV storms from ANN 2-3-1 Model During 

Training 



Figure B-10: Scatter Plot for Category-IV storms from ANN 2-3-1 Model During 

Testing 
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Predicted Streamflow (cumec cp‘ Predicted Streamflow (cumec 
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re B-11: Scatter Plot for Category-IV storms from ANN 2-4-1 Model During 

Training 
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Figure B-12: Scatter Plot for Category-IV storms from ANN 2-4-1 Model During 

Testing 
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Predicted Streamflow (cumec) (w* Predicted Streamflow (cumec) 
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ire B-13; Scatter Plot for Category-IV storms from ANN 3-3-1 Model During 

Training 
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Figure B-14: Scatter Plot for Category-IV storms from ANN 3-3-1 Model During 

Testing 
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Observed Streamflow (cumec) 

Figure B-15: Scatter Plot for Category-IV storms from ANN 3-10-1 Model During 

Training 
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Figure B-16: Scatter Plot for Category-IV storms from ANN 3-10-1 Model During 

Testing 



115 



-4m 




Figure B-18-2: Plot of Goal and Difference Between Training and 
Testing AARE (training S6) from ANN 3-4-1 Model 
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Figure B-19-2: Plot of Goal and Difference Between Training and 
Testing AARE (training S6) from ANN 2-3-1 Model 
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Table B-1; Average Rainfall, Streamflow and Effective rainfall from green-Ampt Method from all 13 Storms 


1 i 

ww 


ON 

cn 


o o o o o o o 


o d 

iS '3 


OOOO^^OOOOOO 


s . 

2 o 
is c: 

1/3 


lo 

ON 

00 

00 

in 

NO 

in 

q 

00 


in 

q 


i> 

r-J 

00 

d 

r-^ 

in 

NO 

cn 

cn 

04 

cn 

Tf 

00 

in 

cn 


W) aj 
2*2 




NO 



NO 

0- 

O' 

Ovl 

rt 

04 

f H 

NO 

04 


o\ 

04 

04 

04 



o o o 


00 CN 


o o o o 


o d 

ill 


o o o o o o o 


a 

d 


in 

CN 

00 

m 

CX3 

04 

in 

in 

00 

2i 

in 

04 

NO 

s 

(U 

o 

cn 

ON 

cn 


cn 

d 

d 


in 

l> 

CM 


C 

r~i 

O 


cn 


cn 

cn 

00 

in 

cn 

04 

cn 







CM 

1— < 


T-l 

r-< 


4) s 
WD - 

2 

<D 


cn 

cs 


o o o o 


o o 


O -3 


u d 

sa 1 

o 

O 

O 

NO 

cn 

o 

o 

o 

o 

O 

O 

o 

1 ^ 


in 

NO 


t". 

1042.2 

ON 

CM 


O'. 

NO 

XT 

r'; 

S o 
is a 
cn 

CM 

CM 

Tt 

04 

CM* 

t> 

CM 

cn 

oo 

NO 

cn 

ON 

CM 

O 

in 

NO 

CM 

od 

04 

d 

CM 

m* 

8 


O 

W) c« 

2‘a 

^ 'd 


NO vq 
O rn 


<N 


Tt 

d 


o o o 


o o 


^ d3 

^ .9 


o o o o o o o 






o o 
is C 

CO 


wn 

CM 

t> 


NO 




Tt" 

cn 


<N 

On 


00 <N 
o 


04 

00 


o d3 
W3 OS 

■<« 


04 oq 
CO d 


o^ 

04 Tt O 
rj- 


NO 

d 


o o o o o 


u d 


I ^ 

b e 

CO 


o 

o 

00 

00 

O' 

00 

ON 

o 

o 

o 



s 

d 

d 




in 

q 

q 

'it 

00 

CM 

q 

cn 

q 

in 

’t 

CM 

't 

in 

ON 

oi 

c^ 

cn 

CM 

5 

CM 

in 


cn 

O' 

CM 


o o o 


0 a 
2| 

1 


^ ^ » 
® ° 8 d 


00 »n 
o| 04 
cn oi 
cn 00 


^ o o o o o o 


M «=! 


o o o 


o o o o o o o 


CO 


II 

CO 


't 

q 

00 

On 

t 3" 

cn 


00 

CM. 

NO 

m 

m 

q 

NO 

NO 

1-H 

in 

NO 

Tj- 

in 

NO 

5 

cn 

cn 

CM 

a\ 

NO 

in 

CM 

NO 

ON 


Tt- 

in 


So 
o 

NO 


® 

WD d 

► 


Csl 

On 


• 0-^000 
o o o 


a 

d 

S a 


»-H 04 cn 


00 On 



stream 

flow 

31.05 

26.676 

22.518 

25.92 

96.66 

124.74 

65.34 

40.5 

33.48 

31.86 

28.89 


Average 

Rainfall 

0 

2.6 

0.4 

33.8 

0 

0 

0 

0 

1.4 

0 

0.4 


Eiiecuve 

Rainfall 

o o o 

1 12.85 

0 

0 

0 

0 

0 

0 

0 

ZIS 

Streamf 

low 



89.64 

89.64 

84.24 

1 105.57 
i 248.94 

380.7 

278.1 

210.6 

164.16 

137.16 

112.32 


^ a 

2 “S 

0) td 
> d 

< ^ 

o o o 

21.6 

0 

0 

0 

0 

0 

0 

2.4 


•^1 
« CS 
04 ‘ds 

SB S 

w cx 

o o o 

10.82 

0 

0 

0 

0 

0 

0 

0 

Sll 

Streamf 

low 

74.52 

66.42 

60.48 

90.72 

1 222.75 

326.7 

226.8 

170.64 

147.15 

128.25 

101.52 


^ a 

2 *S 
S *1 
< 

o o o 

31.6 

0 

0 

0 

2.2 

2.2 

0 

1 


«2 

a.s 

Sd « 
w ^ 

o o o 

0 

37.19 

0 

0 

0 

0 

0 

0 

SIO 

ii 

cn 

202.5 

199.8 

214.65 

272.7 

837 

1042.2 

936.9 

502.2 

267.84 

202,5 

162 


Average 

Rainfall 

0 

0.6 

3.6 

0 

42 

0.4 

0 

0 

0 

1.8 

0 


Effective 

Rainfall 

1 

1 Q 

0 

20.96 

0 

0 

0 

0 

0 

0 

1.76 

0 

ON 

cn 

Streamf! 

ow 

104.76 
; 96.66 
115.02 

442.8 

626.4 

380.7 

294.3 

234.9 

217.08 

191.7 

162 


a 

2*2 
0) .d 

> d 

-<3 P6i 

0 

, 0.6 

27.6 

1.6 

0 

0 

0 ' 
0.4 

1.2 

8.4 

0 


0) 

U'S 

SO! 

o o c 

20.04 

0 

0 

0 

0 

0 

0 

0 

oc 

cn 

Streamf 

low 

248.94 

222.75 

1 906 55 

299.7 

464.4 

572.4 

448.2 

394.2 

337.5 

283.5 

237.6 


« 5- 

2*£ 

‘5 

0.4 

0 

0.2 

28.6 

0 

4.2 

0 

0.2 

0.4 

0 

0 


04 

tSS 

g.E 

! o o c 

0 

0 

0 

0 

0 

0 

0 

9*91 

o 

c/ 

I 

b c 

cn 

221.13 

196.56 

202.5 

426.6 

502.2 

399.6 

340.2 

288.9 

251.1 

224.37 

208.98 


8d^ 

j,.. 

-«:p 

|®dc4JQt'i®®°.-5d° 


« 'c 

^1 



J ! 

fiH 'T 

y|r-^cMcn't*nNot'-ooON2Zl 



References 


Anderson, H. G., and Burt, T. P. (1985). “Hydrological forecasting”. John Willey & Sons 
Publications. 

Anmala, J., Zhang, B., and Govindaraju, R,, S. (2000) “Comparison of ANNs and 
Empirical Approaches for Predicting Watershed Runoff “Journal of Water 
Resources Planning and Management, Vol. 126, No. 3, 156-166 

ASCE task committee (2000). “ Application of neural network in hydrology- 1&2” 

Basha, H. A. (2000). “ Simple nonlinear rainfall-runoff model”. Journal of hydrologic 
engineering ASCE 5( 1 ):25-32. 

Bergstrom, S., and Graham, L.P. (1998). “On the scale problem in hydrological 
modeling”. Journal of Hydrology 211:253-265. 

Beven, K.J. (1989). “Changing ideas in hydrology - the case of physically based 
models”. Journal of Hydrology 105: 157-172. 

Beven, K.J., and Kirkby, M.J. (1979). “A physically based, variable contributing area 
model of basin hydrology”. Hydrological Sciences Journal 24: 43-69 

Binley, A., Beven, K.J., and Elgy, J. (1989). “A physically based model of heterogeneous 
hillslopes. Effective hydraulic conductivities”. Water Resources Research 25(6): 
1219-1226 

Bishop, C. M. (1994). “Neural networks and their applications.’’ Rev. Sci.Instrum., 65, 
1803-1832. 

Bloschl, G. and Sivapalan, M., (1995). “Scale issues in hydrological modelling: a 
review”. Hydrological Processes 9: 251-290 

Buchtele, J., Elias, V., Tesar, M. and Herrmann, A., (1996). “Runoff components 

simulated by rainfall-runoff models”. Hydrological Sciences Journal 41: 49-60 

Bumash, R. J. C., Ferral, R. L. and McGuire, R. A. (1973). “A generalized streamflow 
simulation system conceptual modeling for digital computers.” National Weather 
Service, California Department of Water Resources. 

Chakraborty, K., Mehrotra, K., Mohan, C. K., and Ranka, S. (1992). “Forecasting the 
behaviour of the multivariate time series using neural networks.’’ Neural 
Networks, 5, 961-970. 

Chiu, C. L., and Huang, J. T. (1970). “Nonlinear time varying model of rainfall-runoff 
relation.” Water Resources Research ,6, 1277-1286. 

Chow, V. T., Maidment. D. R., and Mays, L. W. (1988). “Applied hydrology”. McGraw- 
Hill, New York. 

Clark, C. O. (1945). “Storage and the Unit Hydxogxsphl'Transactions of the American 
Society of Civil Engineers, Vol 1 10, pp 1419-1446 

Clarke, R. T. (1973). “A review of some mathematical models used in hydrology, with 
observations on their calibration and use.” Journal of Hydrology, 19, 1-20. 

. Crawford, Norman H., and Linsley, Ray K. (1966). “Digital Simulation in Hydrology: 
Stanford Watershed Model IV,” Civil Engineering Technical Report No. 39, 
Stanford University, Palo Alto, CA. 

Dandy, G. and Maier, H. (1993). “Use of artificial neural networks for forecasting water 


120 



quality in stochastic and statistical methods in hydrology and environmental 
engineering”. An international conference in honor of Professor T. E. Unny, 
University of Waterloo, Ontario. Canada, June. 

Dawdy, D. R. (1982). “A review of deterministic surface water routing in rainfall-runoff 
models. Rainfall-runoff relationship, V. P. Singh, ed.” Water Resource 
Publications, Littleton, Colo., 23-36. 

Dibike Y.B., Solomatine D., Abbott M.B. (1999). “On the encapsulation of numerical 
hydraulic models in artificial neural network”. Journal of Hydraulic Research, 
Vol. 37 (2)- pp. 147-161. 

Dolling, O.R., and Varas, E.A. (2001). “Artificial neural network for Streamflow 
prediction”. Journal of Hydraulic Research, Vol. 40(5), 547-553. 

El-Jabi, N., and Sarraf, S. (1991). “Effect of maximum rainfall position on rainfall-runoff 
Relationship”. Journal of Hydraulic Engineering, 1 17(5), 681-685. 

Elshorbagy, A. and Simonovic, S.P., 2000, Performance evaluation of artificial neural 
networks for runoff prediction, J. Hydrologic Engineering, 5(4), 424-427. 

French, M. N., Krajewski, W. F., and Cuykendall, R. R. (1992). "Rainfall forecasting in 
space and time using a neural network." Journal of Hydrology, 137, 1-31. 

Govindaraju, R. S. (2000). “Artificial neural network in hydrology”. Journal of 
Hydrologic Engineering. , ASCE, 5(2), 115-136. 

Grayson, R. B., Moore, I. D., and McMahon, T. A. (1992). “Physically based hydrologic 
modeling — 2. Is the concept realistic?” Water Resources Research, 28(10), 
2659-2666. 

Gupta, V. K., and Sorooshian, S. (1985). “The relationship between data and the 
precision of parameter estimates of hydrologic models.” Journal of Hydrology., 
81,57-77. 

HEC-1, (late 1960s). U. S. Army Corps of Engineers. Intemate page 

Hjelmfelt, T. A. and Wang, M. (1993). “Artificial neural networks as unit hydro graph 
applications, in C. Y. Kuo (ed.)” Proceedings of the symposium on Engineering 
Hydrology, ASCE, San Francisco, 51 7-522. 

Hjelmfelt, T. A. and Wang, M. (1996). “Predicting Runoff Using Artificial Neural 
Networks,” Surface Water Hydrology, 233-244. 

Homberger, G.M. and E.W. Boyer. (1995). “Recent advances in watershed modeling”. 
U.S. National Report to lUGG, 1991-1994, Reviews of Geophysics, Supplement,' 
pp 949-957 

Horton, R. E. (1940). “An approach towards a physical interpretation of infiltration 
Capacity”. Soil Sci. Soc. Am. Proc., 5, 399-417. 

Houghton-Carr, H. A. (1999) “Assessment criteria for simple conceptual daily rainfall- 
runoff models”. Hydrologic Science Journal 44, 237-261. 

Hsu K.L.,Gupta H.V., Sorooshian S. (1993). “Artificial Neural Network Modeling of the 
Rainfall-Runoff Process”. Water Resources Research Vol. 29(4)pp. 1185-1194. 

Hsu, K., Gupta, V.H. and Sorooshian, S. (1995). Artificial neural network modeling of 
the rainfall-runoff process. Water Resources Research, 31(10), 2571-2530. 

Hsu K, H. V. Gupta, and S. Sorooshian, (1998). “ Streamflow Forecasting Using 
Artificial Neural Networks,” ASCE Water Resources Engineering Conference 
’98, 967-972. 

Jain, Ashu and Indurthy, S. K. V. (2003). “Comparative analysis of event-based rainfall 


121 



runoff modeling techniques-deterministic, statistical and artificial neural 
networks”. Journal of Hydrologic Engineering, 8(2), 93-98. 

Jain, A., and Ormsbee, L. E. (2002). “Evaluation of short-term water demand forecast 
modeling techniques: Conventional v/s artificial intelligence.”/. Am. Water 
Works Assoc., 94(7), 64-72. 

Jain, A., Varshney, A. K., and Joshi, U. C. (2001). “Short-term water demand forecast 
modeling at IIT Kanpur using artificial neural networks.” Water Resour. 
Manage., 15(5), 299-321. 

Jain, S.K. and Chalisgaonkar, D., 2000, Setting up stage-discharge relations using ANN, 
Journal of hydrologic engineering, 5(4), 428-433. 

Jakeman, A. J., and Homberger, G. M. (1993). “How much complexity is warranted in a 
rainfall-runoff model?” Water Resources Research, 29,2637- 2649. 

Kang, K. W., Kim, J. H., Park, C. Y., and Ham, K. J. (1993). “Evaluation of hydrological 
forecasting system based on neural network model.” Proc., 25th Congress oflnt. 
Assoc, for Hydr. Res., International Association for Hydraulic Research, Delft, 
The Netherlands, 257-264. 

Karunanithi, N., Grenney, W.J., Whitley, D., and Bovee, K., 1994, Neural networks for 
river flow predictions. Journal of Computing in Civil Engineering, 8(2), 201-218. 

IChan, H., “Conceptual Modelling of rainfall-runoff systems”, M. Eng. Thesis, National 
University of Ireland, Galway, 1986. 

Kiitanidis, P.K., and R.L. Bras, (1980a-b). “Real-Time Forecasting with a Conceptual 
Hydrologic Model, Analysis of Uncertainty,” Water Resources Research, 16(6), 
1025-1044. 

Kleme, V., (1986). “Operational testing of hydrological simulation models”. 

Hydrological Sciences Journal 31: 13-24 

Konikow, L.F. and Bredehoeff, J.D., (1992). “Ground-water models cannot be 
validated”. Advances in Water Resources 15: 75-83 

t 

Kuczera, G., (1983). “Improved parameter inference in catchment models, 1. Evaluating 
parameter uncertainty”. Water Resources Research 19(5): 1151-1162 

Kull, D., W. and Feldman, A., D. (1998). “Evolution of Clark’s Unit Graph Method to 
Spatially Distributed Runoff “ Journal of Hydrologic Engineering, Vol. 3, No. 1, 
pp. 9-19 

Linsley, R. K. (1982). “Rainfall-runoff models-an overview. Rainfall- Runoff 

relationship, V. P. Singh, ed.”. Water Resources Publications, Littleton, Colo., 3 

22 . 

Linsley, R.K., (1986). “Flood estimates: how good are they?”. Water Resources Research 
22(9): 159S-164S 

League, K.M. (1982). “A comparison of techniques used in rainfall-runoff modeling. 
M.Sc. thesis. 258 pgs. ”, University of British Columbia, Vancouver. 

Markus, M., Salas, J. D., and Shin, H. (1995). “Predicting streamflow based on neural 
Networks”. ASCE, First Inter. Conf. on Water Resources Engineering, San 
Antonio, Texas. 

Martinez, W. L., and Martinez, A. R. (2002). “Computational statistics handbook with 
MATLAB”, Chapman & Hall, London. 


122 



McCulloch, W. S. and Pitts, W. H.- (1943). “A logical calculus of the ideas immanent in 
neural nets”. Bull. Math. Biophys., 5, 115-133. 

Minns, A. W., and Hall, M. J. (1996). “Artificial neural networks as rainfall 
modth'\Hydrologic Science Journal, 41(3), 399-417. 

Mizumura, K., and Chiu, C. L., (1985). Prediction of combined snow-melt and rainfall 
runoff. Journal of Hydrologic Engineering , ASCE, 1 1 1 (2), 179- 193 
Mizumura, K. (1995). “Runoff prediction by simple Tank model, using recession curves” 
Journal of Hydrologic Engineering, ASCE, Nov 1995, 812-818. 

Mohseni, O. and Stefan, H.G., (1998). “A monthly streamflow model”. Water Resources 
Research 34(5): 1287-1298 

Moore, R.D., (1993). “Application of a conceptual streamflow model in a glacierized 
drainage basin”. Journal of Hydrology 150:151-168 
Nandakumar, N. and Mein, R.G., (1997). “Uncertainty in rainfall-runoff model 

simulations and the implications for predicting the hydrological effects of land- 
use change”. Journal of Hydrology 192: 211-232 
Nash, J.E. and B.I. Barsi, “A hybrid model for flow forecasting on large catchments”. 
Journal of Hydrology, 65: 125-137, 1983. 

Nash, J. E. and Sutcliffe, J. V. (1970). “River flow forecasting through conceptual 
models. Part 1, a discussion of principles”. Journal of Hydrology, 10(3), 282-290. 
Nash, J.E. and J.J. Foley, “Linear models of rainfall-runoff systems. In: ‘Rainfall-Runoff 
Relationship’, Proceedings of the International Symposium on Rainfall-Runoff 
modeling. Mississippi State University, May 1981, USA, Edited by V.P. Singh”, 
Water Resources Publications pp.: 51-66, 1982. 

O’Connell, P.E., J.E. Nash, and J.P. Farrell, “River Flow forecasting through conceptual 
models. Part 2. The Brosna catchment at Ferbane”, Journal of Hydrology., 10: 
317-329, 1970. 

O'Connell, P.E. and Todini, E., (1996). “Modeling of rainfall, flow and mass transport in 
hydrological systems: an overview”. Journal of Hydrology 175: 3-16 
O'Loughlin, A., Huber, W., and Chocat, B. (1996). “Rainfall-runoff processes and 
Modeling”. J. Hydr. Res., Delft, The Netherlands, 34, 733- 751. 

Refsgaard, J.C., (1996). “Terminology, modelling protocol and classification of 

hydrological model codes (chapter 2). In; M.B. Abbott and J.C. Refsgaard (eds.). 
Distributed hydrological modelling”. Water Science and Technology Library, 
Vol.22, Kluwer Academic Publishers, Dordrecht, The Netherlands, 321 pp., 17-39 
Refsgaard, J.C., (1997). “Parameterisation, calibration and validation of distributed 
hydrological models”. Journal of Hydrology 198: 69-97 
Rudra P. (1998). “Getting Started with MATLAB 6: A Quick Introduction for Scientists 
and Engineers”. Chapman & Hall, London 
Seibert, J., (1999). “Regionalisation of parameters for a conceptual rainfall-runoff 

model”. Agricultural and Forest Meteorology. 

Shamseldin, A. Y. (1997). “Application of a neural network technique to rainfall-runoff 
modeling”. Journal of Hydrology, 199, 272- 294, 1997. 

Sherman, L. K. (1932). “Streamflow from rainfall by a unit hydrograph method.’’ Eng. 
News Rec., 108, 501-505. 

Singh, V.P. (ed.), (1995). “Computer models of watershed hydrology”. Water Resources 


123 



Publications, Highlands Ranch, Colorado, U.S.A., 1130 pp. 

Singh, V. P. (1988). “Hydrologic systems: Rainfall-runoff modeling”, Vol. I, Prentice- 
Hall, Englewood Cliffs, N.J. 

Singh, S., K. (2004). “Simplified Use of Gamma-Distribution/Nash Model for Runoff 
Modeling” Journal of Hydrologic Engineering, Vol. 9, No. 3, pp. 240-243 

Smith, J., and Eli, R. N. (1995). “Neural network models of rainfall-runoff process”. 

Water Resour. Ping. And Mgmt., ASCE, 121(6), 499-508. 

Sorooshian, S., Gupta, V. K. and Fulton, J. L. (1983). “Evaluation of maximum 
likelihood estimation techniques for conceptual rainfall-runoff models: influence 
of calibration data variability and length on model credibility”. Water Resources 
Research, 19(1), 251-259. 

Subramanaya, K. “Engineering Hydrology”. McGraw-Hill, New York 

Tingsanchali, T., Gautam, M.R., 2000, Application of tank, NAM, ARMA and neural 
network models to flood forecasting. Hydrological Processes, 14(14), 2473-2487 

Todini, E. (1988b). “Rainfall runoff modelling: Past, present and future.” Journal of 
Hydrology, 100, 341-352. 

Todini, E. (1996). “The ARNO rainfall-runoff model”. Journal of Hydrology, 175, 
339-382 

Tokar, A. S., and Johnson, P. A. (1999). “Rainfall-runoff modeling using artificial neural 
networks. ASCE, Journal of Hydrologic Engineering, 4(3), 232-239. 

Tsang, C. F., (1991). “The modeling process and model validation”. Ground Water 29(6): 
825-831 

Werbos, P. (1974). “Beyond regression: new tools for prediction and analysis in the 
behavioral sciences,” PhD dissertation. Harvard University, Cambridge, Mass. 

Xu, C. Y, and Vandewiele, G.L., (1995). “Parsimonious monthly rainfall-runoff models 
for humid basins with different input requirements”. Advances in Water 
Resources 18: 39-48 

Yang, X. and Michel, C. (2000). “Flood forecasting with a watershed model: a new 
method of parameter updating”. Hydrological Sciences Journal, 45(4), 537-546. 

Yapo, P. O., Gupta, V. H., and Sorooshian, S. (1996). “Automatic calibration of 
conceptual rainfall-runoff models: sensitivity to calibration data.” Journal of 
Hydrology, 181, 23-^8. 

Ye, W., Bates, B. C., Viney, N. R., Silvapan, M. & Jakeman, A. J. (1997). “Performance 
of conceptual rainfall-runoff models in low-yielding ephemeral catchments”. 
Water Resources Research, 33, 153-166. 

Zhang B. and R. S. Govindaraju, (1998), “Using Modular Neural Networks to Predict 
Watershed Runoff,” ASCE Water Resources Engineering Conference’98, 897- 
902. 

Zhang, B., and Govindaraju, S. (2000). “Prediction of watershed runoff using Bayesian 
concepts and modular neural networks.” Water Resources Research, 36(3), 753- 
762. 

Zhu, M. Fujita M. (1993) “Application of neural networks to runoff forecast, in K. W. 
Hipel (eds.), stochastic and statistical methods in hydrology and environmental 
Engineering”. An international conference in honor of Professor T. E. Unny, 
University of Waterloo, Ontario. Canada, 3. 


124 



